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PREFACE 


This  is  the  final  technical  report  summarizing  the  adaptive-wall 
wind-tunnel  research  carried  out  in  the  Aerodynamic  Research  Department  of 
the  Calspan  Advanced  Technology  Center,  Buffalo,  NY,  under  Contract  Number 
N00014-77-C-0052,  Task  No.  NR  061-199,  during  the  performance  period  from 
1  November  1976  to  31  October  1980.  The  research  was  sponsored  jointly  by 
the  Office  of  Naval  Research,  with  Mr.  Morton  Cooper  as  technical  monitor, 
and  by  the  Air  Force  Office  of  Scientific  Research,  initially  with 
Mr.  Milton  Rogers  and  later  with  Dr.  James  D.  Wilson  as  technical  monitors. 


The  late  Mr.  R.  J.  Vidal  was  principal  investigator  until  May  1978 
and  was  followed  by  Dr.  J.  C.  Erickson,  Jr.  Dr.  A.  Ritter,  Head  of  the 
Aerodynamic  Research  Department,  had  overall  cognizance  of  the  research. 

The  authors  of  this  final  report  are  indebted  to  several  other  present 
and  former  Calspan  ATC  personnel  who  participated  in  various  tasks  of  this 
research,  namely  P.  A.  Gatlin,  J.  T.  Curtis,  D.  C.  Daughtry,  J.  Nemeth,  Jr., 
J.  P.  Nenni,  R.  E.  Phibbs  and  the  late  A.  F.  Gretch,  The  authors  also  wish 
to  acknowledge  the  valuable  contributions  of  Prof.  W.  R.  Sears  of  the 
University  of  Arizona,  who  was  a  consultant  to  Calspan  ATC  on  this  research 
effort . 
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NOMENCLATURE 


u,  V 


u[v],  v[u] 


x 


0 L 

A 


pressure  coefficient  on  airfoil  surface,  (p-p^  )/<{M 
airfoil  chord 

distance  from  tunnel  centerline  to  control  surface 
Mach  number 
static  pressure 
dynamic  pressure 

Reynolds  number  based  on  chord  length,  c 
free-stream  velocity 

perturhation  velocity  components  in  streamwise, 
normal  directions,  respectively 

u  and  v  found  from  the  exterior-flow  functional - 
relationship  evaluations,  using  v  and  u,  respectively, 
as  boundary  conditions 

streamwise  coordinate  with  origin  at  the  junction 
between  plenum  chambers  6  and  7,  and  16  and  17, 
see  Fig.  1 

angle  of  attack 


Subscript 

free-stream  conditions 
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Section  1 
INTRODUCTION 

The  severity  of  transonic  wind-tunnel  wall  interference  and  the 

inability  of  conventional  procedures  to  correct  the  data  because  of 

the  inherent  nonlinearity  of  the  flow  again  raised  concern  about  ten  years 

ago.  At  that  time,  there  was  a  renewed  importance  attached  to  achievement 

of  reliable  experimental  data  on  flight  vehicle  configurations,  especially 

in  transonic  flow.  Recognition  of  the  deficiencies  of  existing  techniques 

1  2 

and  procedures  led  Ferri  and  Baronti  and  Sears  ,  independently,  to  propose 
the  concept  of  an  adaptive-wall  wind  tunnel  in  which  interference  could  be 
reduced  greatly  at  least,  if  not  eliminated  for  practical  purposes. 

The  reduction  of  interference  in  an  adaptive-wall  wind  tunnel  is 
achieved  by  controlling  the  flow  field  in  the  vicinity  of  the  tunnel  walls. 
Measurement  is  made  of  the  components  of  the  disturbance  velocity  at  discrete 
points  along  imaginary  control  surfaces,  or  interfaces,  in  the  flow  field 
within  the  tunnel.  A  theoretical  representation  for  the  flow  field  external 
to  the  control  surfaces,  including  the  boundary  condition  for  unconfined 
flow,  i.e.,  that  all  disturbances  vanish  at  infinity,  is  used  to  determine 
if  those  measured  velocity  components  satisfy  functional  relationships  which 
are  consistent  with  interference-free  flow.  If  they  do  not,  an  iteration 
procedure  provides  a  new  approximation  for  the  flow  field  at  the  interfaces, 
and  the  flow  control  in  the  vicinity  of  the  walls  is  readjusted  successively 
until  the  measured  quantities  are  consistent  with  the  boundary  condition  for 
unconfined  flow.  In  this  way,  theory  and  experiment  are  combined  to  minimize 
wall  interference.  The  concept  of  an  adaptive  wall-tunnel  is  discussed  more 
completely  in  References  2  to  5. 

A  program  of  research  has  been  in  progress  at  Calspan  since  1971 
to  develop  and  demonstrate  the  feasibility  of  the  adaptive-wall  concept.  The 
concept  was  demonstrated  theoretically  in  1973  and  1974  by  iterating  idealized 

*  Ferri,  A.  and  Baronti,  P.  "A  Method  for  Transonic  Wind  Tunnel  Corrections" 
AIAA  Journal,  Vol.  11,  No.  1,  pp.  63-66,  January  1973. 

2 

Sears,  W.R,  "Self-Correcting  Wind  Tunnels"  (The  Sixteenth  Lanchester 
Memorial  Lecture)  Calspan  Report  No.  RK-5070-A-2,  July  1973;  also  the 
Aeronautical  Journal,  Vol.  78,  No.  7S8/759,  pp.  80-89,  February/March  1974. 
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numerical  simulations  of  the  flow  in  two-dimensional  adaptive-wall  tunnels 

3  4 

for  incompressible  and  transonic,  supercritical-wall  conditions  .  Also 

in  1973,  a  two-dimensional  adaptive-wall,  transonic  test  section  was 

designed,  fabricated  and  installed  in  the  circuit  of  the  Calspan  One-Foot 

Tunnel.  The  test  section,  which  is  described  more  fully  in  Section  3 

and  in  Refs.  4  and  5,  consists  of  perforated  upper  and  lower  walls  with 

segmented  plenum  chambers.  Each  plenum  has  an  individual  pressure  control 

to  provide  either  suction  or  blowing.  A  model  with  an  NACA  0012  airfoiL 

4  6 

section  and  a  6- inch  chord  was  fabricated  and  tested  *  in  the  Calspan 
Eight-Foot  Transonic  Wind  Tunnel  to  establish  the  airfoil  characteristics 
in  unconfined  flow  at  a  Reynolds  number,  Re  of  1.00  x  106  based  on  chord 

C  i 

length,  c.  Aerodynamic  data  were  obtained  at  free-stream  Mach  numbers, M  , 
from  0.40  to  0.95  and  angles-of-attack,  ,  from  -2°  to  8®. 

4  5 

An  experimental  program  ’  was  carried  out  during  1974  and  1975 
in  the  One-Foot  Tunnel  for  flows  which  were  supercritical  at  the  model, 
but  subcritical  at  the  control  surfaces  and  walls,  i.e.  for  M  0.725. 

Initial  experiments  with  the  tunnel  operated  so  as  to  simulate  conventional 

4 

perforated-wall  tunnel  operation  displayed  significant  wall  interference 

for  these  free-stream  conditions,  since  the  solid  blockage  of  the  6-inch 

chord  model  is  6%.  Iteration  at  several  flow  conditions  in  this  regime  of 
4  5 

operation  indicated  ’  satisfactory  convergence  to  unconfined  flow  and 
elimination  of  the  wall  interference. 

3 

Erickson,  J.C.,  Jr.  and  Nennl,  J.P.  "A  Numerical  Demonstration  of  the 
Establishment  of  Unconfined-Plow  Conditions  in  a  Self-Correcting  Wind  Tunnel" 
Calspan  Report  RK-5070-A-1,  November  1973. 

4 

Vidal,  R.J.,  Erickson,  J.C.,  Jr.  and  Catlin,  P.A.  "Experiments  with  a  Self- 
Correcting  Wind  Tunnel"  AGARD-C'P-174 ,  October  1975;  also  Calspan  Report. 

No.  RK-5070-A-4 ,  October  1975. 

®  Sears,  W.R.,  Vidal,  R.J.,  Erickson,  J.C.,  Jr.  and  Ritter,  A.  "Interference- 
Free  Wind-Tunnel  Flows  by  Adaptive-Wall  Technology"  ICAS  Paper  No.  76-02, 
10th  Congress  of  the  International  Council  of  the  Aeronautical  Sciences, 
Ottawa,  Canada,  3-8  October  1976;  also  Journal  of  Aircraft,  Vol .  14,  No.  11, 
pp.  1042-1050,  November  1977. 

6  Vidal,  R.J.,  Catlin,  P.A.  and  Chudyk,  D.W.  "Two-Dimensional  Subsonic 
Experiments  with  an  NACA  0012  Airfoil"  Calspan  Report  No.  RK-5070-A-3, 
December  1973. 
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The  success  achieved  at  these  subcritical-wall  conditions  led  . 
to  extension  of  the  research  effort  to  encompass  tunnel  operating  conditions 
at  higher  Mach  numbers  for  which  shock  waves  generated  by  the  model  reach 
the  control  surfaces  and  walls.  Moreover,  consideration  was  given  to  a 
perforated- wall ,  segmented-plenum  implementation  of  the  adaptive-wall 
concept  in  three  dimensions.  This  report  summarizes  the  research  performed 
to  meet  these  objectives  with  joint  sponsorship  by  the  Office  of  Naval 
Research  (ONR)  and  the  Air  Force  Office  of  Scientific  Research  (AFOSR) . 

There  have  been  three  principal  tasks  in  the  investigation. 

First,  experimental  research  was  carried  out  in  the  two-dimensional  test 

section  of  the  One-Foot  Tunnel.  During  1976  and  1977,  this  task  complemented 
7  8 

and  supplemented  the  experiments  performed  with  sponsorship  by  the 
Arnold  Engineering  Development  Center  (AEDC) .  During  1980,  additional 

9 

experiments  were  carried  out  as  a  continuation  of  research  accomplished 
during  1978  and  1979  with  further  AEDC  sponsorship.  The  second  task, 
performed  during  1978  and  1979,  was  the  development  of  more  realistic 
numerical  simulations  of  an  airfoil  in  a  two-dimensional  adaptive-wall  test 
section,  including  the  influence  of  tunnel-wall  boundary  layers.  The  third 
task  was  a  preliminary  investigation  of  the  design  requirements  for  three- 
dimensional  adaptive-wall  test  sections,  with  particular  attention  given  to 
the  adaptability  of  the  Calspan  One-Foot  Tunnel  to  meet  these  requirements. 
This  task  was  performed  during  1978  and  1979. 

7  Vidal,  R.J.  and  Erickson,  J.C.,  Jr.  "Experiments  on  Supercritical  Flows 
in  a  Self-Correcting  Wind  Tunnel"  AIAA  Paper  No.  78-788,  AIAA  10th 
Aerodynamic  Testing  Conference,  San  Diego,  California,  19-21  April  1978. 

®  Vidal,  R.J.  and  Erickson,  J.C.,  Jr.  "Research  on  Adaptive-Wall  Wind 
Tunnels"  AEDC  Report  No.  AEDC -TR- 78- 36,  November  1978. 

^  Erickson,  J.C.,  Jr.,  Wittliff,  C.F. .  and  Daughtry,  D.C.  "Further 
Investigations  of  Adaptive -Wall  Wind  Tunnels"  AEDC  Report  No. 
ABDC-TR-80-34,  October  1980. 


For  completeness  and  to  place  in  perspective  the  progress  on 
adaptive-wall  wind  tunnels,  a  brief  survey  of  adaptive-wall  experiments 
at  other  laboratories  is  given  in  Section  2.  The  adaptive-wall  test 
section  and  airfoil  models  for  the  Calspan  One-Foot  Tunnel,  including 
modifications  that  were  made  during  the  present  investigation,  are  described 
in  Section  S.  Also  described  in  this  Section  is  the  development  of  instru¬ 
mentation  and  calibration  techniques  for  measuring  the  streamwise  disturbance 
velocity  component,  u  ,  and  the  normal  component,  v.  In  Section  4,  results 
of  experimental  iterations  toward  interference- free  flow  conditions  with 
a  41-blockage  NACA  0012  airfoil  model  are  presented  for  supercritical  flow 
conditions  at  the  walls.  These  experiments,  which  were  performed  during 
1980  at  ■  0.9  and  (X  =  3°,  2°  and  1°,  have  achieved  reasonable  success 
and  have  extended  the  applicability  of  the  adaptive-wall  concept  to  more 
severe  transonic  conditions  than  had  been  accomplished  previously. 

Appendix  I  contains  a  list  of  papers,  reports  and  presentations 
that  were  prepared  during  this  investigation.  Appendix  II  contains  a  summary 
of  the  Calspan  Design  studies  for  a  three-dimensional  adaptive-wall  test 
section.  Finally,  a  copy  is  attached  of  AIAA  Paper  No.  81-0160  entitled, 
Numerical  Simulations  of  a  Segmented-Plenum,  Perforated,  Adaptive-Wall  Wind 
Tunnel .  This  paper  by  J.C.  Erickson,  Jr.  and  G.F.  Homicz  was  presented  at 
the  AIAA  19th  Aerospace  Sciences  Meeting  in  St  Louis  on  12-15  January  1981, 
and  describes  the  results  of  the  work  on  the  second  task. 


Section  2 


SURVEY  OP  ADAPTIVE- WALL  EXPERIMENTAL  RESEARCH 

There  have  been  several  different  approaches  to  practical  wall- 
control  methods  for  adaptive-wall  wind  tunnels.  To  date,  with  exceptions 
to  be  noted,  the  experimental  research  has  been  carried  out  in  two-dimensional 
flows.  Two  fundamentally  different  wall-control  methods  have  been  pursued, 
namely  changing  the  shape  of  impermeable,  flexible  walls  and  using  ventilated 
walls,  either  perforated  or  slotted,  with  segmented  plenum  chambers  and/or 
porosity  control.  In  this  section,  a  brief  survey  is  given  of  the  experimental 
results  that  have  been  achieved  elsewhere.  There  will  be  no  discussion  here 
of  any  theoretical  efforts  in  support  of  the  concept. 

Impermeable,  flexible-wall  tunnels  have  been  developed  by  three 
European  groups,  namely  thosed  led  by  Chevallier  at  ONERA/Chalais ,  Goodyer 
at  the  University  of  Southampton,  and  Ganzer  at  the  Technical  University 
of  Berlin.  In  all  of  these  implementations,  the  wall  positions,  with 
corrections  for  boundary-layer  displacement,  are  used  to  determine  the 
velocity  component,  v  ,  normal  to  each  control  surface  (the  wall  in  this 
case)  and  static  pressure  measurements  along  the  walls  are  used  to  determine 
the  streamwise  component,  u  . 

Chevallier10  achieved  successful  iteration  to  unconfined  flow  in 
1975  with  a  4.4%-blockage  NACA  64A010  model  for  Mw  up  to  0.85.  This  tunnel 
configuration  and  associated  techniques  have  been  adapted  in  a  two- 
dimensional  tost  section  of  0.4  m  x  0.4  m  for  the  T2  Transonic  Tunnel  at 
ONERA/CERT  at  Toulouse. 

10  Chevallier,  J.P.  "Soufflerie  Transsonique  a  Parois-Adaptables"  AGARD- 
CP-174,  October  1975;  also  translated  into  English  as  European  Space  Agency 
Report  ESA-TT-326,  October  1976,  available  as  NASA  Accession  No.  N77-13085. 

11  Poisson-Quinton,  P.  "Some  New  Approaches  for  Wind-Tunnel  Testing  Through 
the  Use  of  Computers"  AIAA  Paper  No.  79-0707,  First  Intersociety  Atlantic 
Aeronautical  Conference,  Williamsburg,  Virginia,  26-28  March  1979;  also 
ONERA  TP  No.  1979-24. 
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Goodyer,  et  al .  have  obtained  similar  success,  initially  at  low 

12  13 

speeds  in  1975,  and  then  at  transonic  speeds  in  1980.  In  these  latter 

experiments,  unconfined- flow  conditions  were  achieved  at  M*,  -  0.89 

and  CL  ■  4°  with  an  8%-blockage  NACA  0012-64  model.  However,  it  should 

be  noted  that  although  shock  waves  extend  to  the  walls,  the  exterior-flow 

calculation  was  based  on  subcritical-flow  theory.  Goodyer' s  configuration 

14 

and  techniques  are  being  implemented  in  a  two-dimensional  test  section  of 
0.33  m  x  0.33  m  in  the  0.3m  Transonic  Cryogenic  Tunnel  at  the  NASA  Langley 
Research  Center. 

Ganzer  achieved  fully-automated  iterative  control  of  the 

walls  for  8%-blockage  models  of  NACA  0012  and  supercritical  CAST  7  airfoils 

up  to  Mm  >0.82  during  1979  and  1980.  Again,  subcritical  exterior-flow 

calculations  were  performed.  Ganzer  also  has  carried  out*^  preliminary 

experiments  in  the  two-dimensional  tunnel  with  a  three-dimensional  model, 
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and  the  other  groups  plan  *  to  do  so  as  well.  The  objective  in  these 
experiments  is  to  achieve  partial  control  and  then  to  use  the  measured 
wall  pressures  to  compute  residual  corrections  along  the  lines  first 

Goodyer,  M.J.  "A  Low-Speed  Self  Streamlining  Wind  Tunnel"  AGARD- 
CP-174,  October  1975. 

Goodyer,  M.J.  and  Wolf,  S.W.D.  "The  Development  of  a  Self-Streamlining 
Flexible  Walled  Transonic  Test  Section"  AIAA  Paper  No.  80-0440, 

AIAA  11th  Aerodynamic  Testing  Conference,  Colorado  Springs,  Colorado, 
18-20  March  1980. 

14 

Ladson,  C.L.  "A  New  Airfoil  Research  Capability"  in  Advanced  Technology 
Airfoil  Research,  Vol.  1,  NASA  CP-2045,  Part  1,  March  1979. 

15 

Ganzer,  U.  "Windkanale  mit  Adaptiven  Wanden  zur  Beseitigung  von 
Wandinterferenzen^'  Zeitschrift  fur  Flugwissenschafteu  und 
We ltraumfors chung,  Vol.  3,  No.  2,  pp.  129-133,  1979;  also  translated  into 
English  as  NASA  TM-75501,  August  1979. 

Ganzer,  U.  "Adaptable  Wind  Tunnel  Walls  for  2-D  and  3-D  Model  Tests" 

ICAS  Paper  No.  23-3,  12th  Congress  of  the  International  Council  of 
the  Aeronautical  Sciences,  Munich,  Germany,  12-17  October  1980. 
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proposed  by  Kemp  '  and  later  treated  by  Capelier,  et  al .  and  Murman 

Ganzer  also  describes1^  plans  for  building  a  three-dimensional  flexible-wall 


test  section. 

Emphasis  in  the  United  States,  namely  at  AEDC  and  NASA  Ames  Research 

Center  besides  Calspan,  has  been  on  ventilated  adaptive-wall  configurations. 

21 

Also,  although  the  complete  adaptive-wall  procedures  were  not  used,  Weeks 
accomplished  some  wall  control  at  the  Air  Force  Flight  Dynamics  Laboratory  in 
1975  by  slotted  walls  with  contoured  slots. 

22 

AEDC  results,  .  published  in  1979,  describe  experiments  with  two 
different  wall  configurations  in  the  Aerodynamic  Wind  Tunnel  (IT)  up  to  Mrt  « 
0.80  for  an  NACA  0012  model  with  6%-blockage.  One  wall  configuration  has  pro¬ 
visions  for  varying  the  hole  angle  of  a  perforated  wall  as  a  function  of  the 
streamwise  direction  by  means  of  multiple  rows  of  bored  spheres  connected  by 
rods  normal  to  the  streamwise  direction.  The  other  wall  configuration  has 
variable-porosity  walls  formed  by  matched  sliding  plates  that  have  inclined 
holes  like  those  in  the  AEDC  Aerodynamic  Wind  Tunnel  (4T) .  These  walls  can  be 
adjusted  globally  to  give  a  porosity  that  is  uniform  on  each  wall,  but  may 
differ  between  the  top  and  bottom  walls.  In  each  implementation,  the  top  and 
bottom  walls  were  both  vented  to  a  single  plenum  chamber,  the  pressure  of  which 


Kemp,  W.B.,  Jr.  "Toward  the  Correctable- Interference  Transonic  Wind 
Tunnel"  AIAA  Paper  No.  76-1794,  AIAA  9th  Aerodynamic  Testing  Conference, 
Arlington,  Texas,  7-9  June  1976. 
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Kemp,  W.B.,  Jr.  "TWINTAN:  A  Program  for  Transonic  Wall  Interference 
Assessment  in  Two-Dimensional  Wind  Tunnels"  NASA  TM-81819,  May  1980. 

19 

Capelier,  C.,  Chevallier,  J-P.  and  Bouniol,  F.  "A  New  Method  for 
Correcting  Wall  Interference"  La  Recherche  Aerospatiale,  1978,  No.  1, 
January-February  1978,  pp.  1-11;  also  translated  into  English  as  European 
Space  Agency  Report  ESA-TT-491,  August  1978,  available  as  NASA  Accession 
No.  N 79-11997. 

20 

Murman,  E.M.  "A  Correction  Method  for  Transonic  Wind  Tunnel  Wall  Inter¬ 
ference"  AIAA  Paper  No.  79-1533,  AIAA  12th  Fluid  and  Plasma  Dynamics 
Conference,  Williamsburg,  Virginia,  24-26  July  1979. 

21 

Weeks,  T.M.  "Reduction  of  Transonic  Slotted-Wall  Interference  by  Means 
of  Slat  Contouring:  Air  Force  AFFDL-TR-74-139,  March  1975. 

22 

Kraft,  E.M.  and  Parker,  R.L.,  Jr.  "Experiments  for  the  Reduction  of 
Wind  Tunnel  Wall  Interferenced  by  Adaptive-Wall  Technology" 

AEDC  Report  No.  AEDC-TR- 79-51,  October  1979. 


was  adjustable  in  magnitude.  In  both  experiments,  u  was  measured  by  static 
pipes  and  v  by  aerodynamic  probes.  Both  subsonic  and  transonic  small- 
disturbance  theory  were  used,  as  appropriate,  to  evaluate  the  functional  rela¬ 
tionships  for  unconfined  flow.  Sufficient  control  existed  with  both  wall 
configurations  to  reduce  the  interference  significantly,  but  complete  iterations 
to  unconfined  flow  were  not  achieved. 

ABDC,  in  a  paper  published  in  1980,  also  obtained  significant  reduc- 
23 

tions  in  interference  in  three-dimensional  flows  for  a  wing/ fuselage/hori¬ 
zontal-tail  configuration  in  Aerodynamic  Wind  Tunnel  (4T) .  In  these  experiments, 
global  porosity  variation  on  the  top,  bottom  and  side  walls  was  used  as  with 
the  two-dimensional  variable-porosity  walls,  again  with  constant  plenum  pressure. 
A  single  translating  probe  was  used  to  measure  both  u  and  v  at  a  control  surface 
of  rectangular  cross  section  that  surrounded  the  model . 

The  NASA  Ames  Rosearch  Center  adaptive-wall  test  section24  has 
slotted  top  and  bottom  walls  with  ten  independently-controlled  plenum 
chambers  behind  each  wall.  Laser  velocimetry  is  used  to  measure  v  at  two 
different  distances  from  the  model  and  the  corresponding  functional- 
relationship  evalutions  were  carried  out  by  subsonic  theory.  Measured 
influence  functions  based  on  plenum  pressure  adjustment  were  used  to 
determine  control  adjustments.  Successful  iteration  to  unconfined  flow  was 
reported  in  1980  up  to  about  Mg,  -  0.8  for  an  NACA  0012  airfoil  model  with 
7%-blockage.  The  use  of  this  two-dimensional  tunnel  for  examining  fully 
three-dimensional  flows  in  the  manner  of  the  European  experiments  also 
has  been  suggested24, 

23 

Parker,  R.L.,  Jr.  and  Sickles,  W.L.  "Application  of  Adaptive  Wall  Techni¬ 
ques  in  a  Three-Dimensional  Wind  Tunnel  with  Variable  Wall  Porosity"  AIAA 
Paper  No.  80-0157,  AIAA  18th  Aerospace  Sciences  Meeting,  Pasadena, 

California,  14-16  January  1980. 

24 

Bodapati,  S.,  Schairer,  E.  and  Davis,  S.  "Adaptive -Wall  Wind-Tunnel 
Development  for  Transonic  Testing"  AIAA  Paper  No.  80-0441,  AIAA  11th 
Aerodynamic  Testing  Conference,  Colorado  Springs,  Colorado,  18-20  March  1980. 


8 


Section  3 


EXPERIMENTAL  FACILITY  AND  INSTRUMENTATION 


3.1  Test  Section  and  Airfoil  Models 

4  5 

The  Calspan  One-Foot  Tunnel  *  ,  shown  schematically  in  Figure  1, 
is  a  closed-circuit,  continuous  flow  facility.  The  two-dimensional  adaptive- 
wall  test  section  is  12-inches  high,  10-inches  wide  and  56-inches  long  with 
perforated  top  and  bottom  walls  of  22.5%  open-area  ratio.  The  plenum 
chambers  behind  the  perforated  walls  have  been  divided  into  18  segments,  10 
on  the  top  and  8  on  the  bottom,  and  each  segment  is  connected  to  a  pressure 
and  a  suction  source  through  individual  control  valves.  Hie  pressure  source 
is  the  tunnel  stilling  chamber,  and  the  suction  source  is  an  auxiliary 
compressor  discharging  into  the  tunnel  circuit  in  the  diffuser.  Six  plenum 
chambers  in  the  immediate  vicinity  of  the  model  have  provisions  for  a 
distributed  porosity  which  can  be  varied  linearly  in  the  streamwise  direction. 
This  capability  was  not  exercised  in  these  experiments  and  a  constant  open- 
area  ratio  was  used. 

As  mentioned  in  the  Introduction,  a  model  with  an  NACA  0012 

airfoil  section  and  a  6-inch  chord  (6%  blockage  in  the  One-Foot  Tunnel) 

4  6 

was  fabricated  and  tested  ’  in  the  Calspan  Eight-Foot  Transonic  Wind  Tunnel. 

This  model  is  instrumented  with  a  row  of  pressure  orifices  on  both  the  upper 

and  lower  surfaces,  and  with  a  three-component  force  balance  supporting  a 

metric  section  on  the  tunnel  centerline.  An  initial  series  of  experiments 

7  8 

with  this  airfoil  model  at  higher  Mach  numbers  had  shown  '  that  the  available 

wall  control  was  limited  so  that  successful  iteration  to  unconfined  flow 

could  not  be  achieved.  This  recurring  inability  to  achieve  full  control  of 

the  flow  field  at  all  plenum  sections  led  to  an  analysis  of  the  flow  in  the 
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auxiliary  suction  system  and  the  test  section  *  . 


Basically,  the  analysis  consists  of  writing  a  pressure  balance 
for  the  auxiliary  suction  system.  The  analysis  includes  the  operating 
characteristics  of  the  auxiliary  compressor,  the  losses  in  the  pipes,  the 
recompression  pressure  loss  arising  from  the  compressor  discharge  into  the 
diffuser,  the  pressure  drop  across  the  perforated  walls,  and  the  required 
unconfined-flow  pressure  at  the  wall  locations  due  to  the  model.  When 
suitable  engineering  approximations  for  these  pressure  terms  were  made, 
using,  in  part,  data  measured  in  the  One-Foot.  Tunnel,  the  predicted  limits 
on  available  control  agreed  reasonably  well  with  experimental  observations. 
One  conclusion  from  this  analysis  was  that  the  recompression  penalty  could 
be  reduced  considerably.  This  was  accomplished  by  introducing  an  area  change 
in  the  tunnel  diffuser  at  the  location  where  the  flow  from  the  auxiliary 
blower  is  vented  into  the  tunnel  circuit. 

The  analysis  also  illustrates  the  three  major  trade-offs  available 
in  the  design  and  application  of  adaptive-wall  test  sections  of  this 
configuration;  namely  compression  ratio,  wall  open-area  ratio  and  model  size. 
Tunnel  performance  could  be  improved  by  increasing  the  compression  ratio, 
but  the  improvement  would  be  available  only  at  Mach  numbers  above  about  0.75, 
because  the  wall  perforations  would  be  choked  at  Mach  numbers  below  that 
value.  That  restriction  could  be  relaxed  considerably  by  using  a  larger 
open-area  ratio,  although  the  shock-wave  reflection  characteristics  might  not 
be  as  favorable.  Decreasing  the  model  size  from  6%  blockage  to  a  smaller 
value  would  decrease  the  magnitude  of  the  disturbance  velocities  at  the 
walls  and  would  improve  the  tunnel  performance.  Consequently,  an  NACA  0012 
airfoil  model  with  a  4-inch  chord  (4%  solid  blockage)  was  constructed.  It 
has  a  row  of  pressure  orifices  along  its  centerline  on  the  upper  and  lower 
surfaces,  but  there  are  no  provisions  for  measuring  the  forces  and  pitching 
moment  directly.  This  model  was  used  in  the  experiments  to  be  described 
below. 


3.2  Initial  Plow-Velocity  Measurement  Techniques 


In  the  experiments  prior  to  those  of  Ref.  9,  local  flow  angle  was 

measured  by  aerodynamic  probes  and  the  local  static  pressure  by  static  pipes. 

These  measurements  were  used,  respectively,  to  determine  v  and  u  .  The 
4  5 

static  pipes,  *  each  with  52  static  orifices  along  its  length,  were  located 
with  their  centerlines  nominally  4.0-inches  above  and  below  the  tunnel 
centerline  and  had  a  0.5-inch  diameter.  Eighteen  flow-angle  probes  ’  were 
mounted  nominally  4.5-inches  above  and  below  the  centerline  with  each  one 
opposite  the  center  of  a  plenum  chamber.  Blockage  and  probe  wake  consider¬ 
ations  limited  the  number  of  probes  to  this  quantity,  which  made  it  difficult 
to  obtain  sufficient  measurements  to  define  adequately  the  v  distributions 
in  the  vicinity  of  the  model. 

Operational  procedures  required  that  the  distributions  of  the 
measured  normal  velocity  component  v  be  used  as  the  boundary  conditions 
on  the  external -flow  calculations  to  provide  the  next  approximation  to  u. 

It  was  difficult  to  carry  out  conventional  interpolation  procedures  accurately 
based  on  the  limited  number  of  v  measurements  available.  Fortunately, 
however,  sufficient  numbers  of  u  measurements  were  available  to  provide 
a  good  definition  of  their  distributions.  Based  on  this  situation,  a 
procedure  was  devised  to  provide  a  better  interpolation  in  the  v  data. 

The  first  step  in  this  interpolation  procedure  was  to  use  the 
extensive  u  data,  as  interpolated  by  a  cubic  spline  (smoothed  or  not,  as 
desired),  to  calculate  the  corresponding  unconfined- flow  distribution  v[u]. 
This  calculation  used  a  finite-difference  solution  of  the  transonic  small- 

g 

disturbance  equations  .  Next,  the  difference  v  -  v[u]  was  determined  at  each 
x  loaation  where  v  measurements  wore  made.  These  differences  were  then  inter¬ 
polated  linearly  in  x,  so  that  the  resultant  interpolation  in  v  was  found  by 
adding  the  interpolated  difference,  v  -  v[u],  to  the  calculated  v[u]  distri¬ 
bution.  Once  the  v  interpolation  had  been  found,  the  exterior-flow  calcula¬ 
tion  of  u[v]  was  carried  out  to  give  the  next  approximation  to  u  which  was 
set  in  the  tunnel.  Again,  transonic  small-disturbance  solutions  were  used. 


This  procedure  thus  used  the  shape  of  the  v[u]  distribution  as 
a  basis  for  interpolation  since  the  distribution  was  well-defined  experi¬ 
mentally.  Clearly, as  convergence  of  the  tunnel  iteration  was  approached, 
the  differences  v  -  v[u]  approached  zero  and  the  interpolation  improved. 
Limited  use  of  this  procedure  in  the  experiments  of  Ref.  8  and  a  few 
subsequent  ones  indeed  showed  an  improvement  in  the  interpolation.  However, 
questions  still  existed  about  the  actual  behavior  of  v  between  measuring 
points.  Moreover,  the  probes  have  limited  the  research  in  other  ways. 

The  small  probes  used  are  very  sensitive  to  contamination  from  oil  present 
in  the  tunnel  air  stream  and  must  be  cleaned  frequently.  This  cleaning  can 
disturb  them  and  render  their  calibrations  for  zero  flow  angle  suspect,  and 
thus  require  frequent  re-calibration.  Therefore,  an  alternative  technique 
was  sought. 

3.3  Revised  Plow-Velocity  Measurement  Techniques 

The  alternative  technique  selected  for  measuring  v  involves 

measuring  the  static  pressures  at  one  control  surface  and  the  difference 

between  those  pressures  and  the  pressures  at  a  second  surface  slightly 

farther  away  from  the  model.  In  effect,  this  can  be  regarded  as  measuring 

the  local  static  pressure  and  its  gradient,  from  which  the  streamwise 

derivative  of  the  normal  velocity,  dv/dx,  can  be  inferred.  The  advantage 

in  this  measurement  technique  is  that  static  pressure  is  easy  to  measure 

with  good  precision  and  one  can  easily  obtain  good  spatial  resolution. 

25 

The  TSFOIL  computer  code  was  used  to  calculate  preliminary  estimates  of 
the  pressure  differences  to  be  expected,  and  these  differences  are  readily 
measurable,  in  the  immediate  vicinity  of  the  airfoil. 

Development  of  a  new  static-pipe  technique  to  achieve  this  was 
performed  with  AEDC  sponsorship  and  is  described  in  Ref.  9.  Briefly,  the 
new  pipes,  of  0.625-inch  diameter,  have  diametrically  opposed  orifices  on 

2S . ' . 

Murman,  E.M.,  Bailey,  F.R.  and  Johnson,  M.L.  "TSFOIL-A  Computer  Code  for 
Two-Dimensional  Transonic  Calculations,  Including  Wind-Tunnel  Wall  Effects 
and  Wave-Drag  Evaluation"  Paper  No.  26  in  Aerodynamic  Analyses  Requiring 
Advanced  Computers,  NASA- SP- 34 7 -Part  2,  March  1975. 
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their  top  and  bottom  in  the  vicinity  of  the  model.  There  are  18  pairs  of 
these  orifices  extending  9  inches  upstream  and  11  inches  downstream  of  the 
junction  between  plenum  chambers  6  and  7,  as  shown  in  Figure  1.  These  dual 
orifices  span  the  region  where  the  static  pressure  differences  are  equal  to, 
or  greater  than,  the  resolution  capability  of  the  pressure  transducers  being 
used.  Upstream  and  downstream  of  this  region,  the  static  pipes  have  orifices 
which  extend  the  full  length  of  the  test  section  along  the  side  of  the  pipe 
facing  the  model.  The  most  forward  static  pressure  orifice  on  each  pipe 
(which  is  located  at  the  beginning  of  the  test  section)  is  connected  to  a 
manifold  and  the  reading  is  taken  to  be  the  free-stream  static  pressure,  p^  . 
All  remaining  33  pressures  on  the  model  sides  of  the  pipes  are  measured 
relative  to  .  In  addition,  the  differential  pressures  between  the  18 
opposing  pairs  of  orifices  on  each  pipe  are  measured.  The  differential 
pressure  transducers  used  have  a  probable  error  of  0.001  psi  or  less,  and 
the  read-out  system  has  a  resolution  of  0.001  psi. 

After  fabrication,  the  new  static  pipes  were  mounted  with  their 
centerlines  nominally  4.0  inches  from  the  test-section  centerline,  as  the 
original  pipes  were.  However,  in  order  to  use  the  probe  flow-angle  data 
in  conjunction  with  the  differential  pressure  data,  the  probes  were  re¬ 
located  to  the  plane  of  the  static  pipe  centerlines  from  their  original 
locations  4.5  inches  from  the  test-section  centerline.  In  the  original 
configuration,  four  of  the  probes  in  the  vicinity  of  the  model  were  mounted 
through  the  test-section  windows  and  their  frames.  It  was  not  feasible  to 
relocate  or  remove  these  probes,  so  they  were  retracted  to  lie  against  the 
side  walls.  However,  two  new  flow-angle  probes  were  fabricated  and  mounted 
in  the  vicinity  of  the  model.  The  resulting  locations  of  the  flow-angle 
probes  at  the  lower  control  surface  are  shown  in  the  side  view  of  Figure  1, 
while  the  upper  probes  are  shown  in  the  top  planform  view.  The  lateral 
staggering  of  the  probes  was  chosen  to  avoid  interference  effects. 
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The  flow  about  static  pipes  in  the  presence  of  mode 1/wall- induced 
9 

disturbances  was  analyzed  within  the  framework  of  conventional  slender-body 

theory  for  subcritical,  compressible  flow  and  the  theory  of  matched  asymptotic 

expansions.  This  analysis  relates  the  average  of  the  pressure  measurements 

across  the  pipe  to  u  at  the  pipe  centerline.  In  addition,  the  analysis 

relates  the  differential  pressure  measurements  across  the  pipe  to  the  stream- 

wise  derivative  of  the  v  at  the  pipe  centerline.  The  derivative,  dv/dx,  can 

be  integrated  to  obtain  the  distribution  of  v  if  at  least  one  independent 

9 

measurement  of  v  is  made.  Based  on  experiences  with  supercritical  flows  in 
which  the  shock  wave  from  the  model  intersects  the  pipe,  separate  integrations 
of  dv/dx  are  carried  out  upstream  and  downstream  of  the  shock,  thus  requiring 
independent  v  measurements  in  each  region. 


As  mentioned  previously,  it  is  necessary  to  calibrate  the  probes 
in  the  tunnel  periodically  to  determine  their  apparent  flow-angle  readings 
for  uniform  flow  parallel  to  the  tunnel  centerline.  The  procedure  for  this 
part  of  the  overall  probe  calibration  technique  is  to  use  the  wall  control 

to  set  a  uniform  static  pressure  along  the  entire  length  of  both  static 
pipes.  The  flow  is  assumed  to  be  parallel  to  the  tunnel  centerline  and  the 
probe  readings  are  taken  to  be  those  for  zero  flow  inclination.  An  efficient 
independent  measurement  technique  to  assure  that  the  flow  is  parallel  to  the 
centerline  is  not  available.  However,  the  procedure  just  described  was 
verified  independently  early  during  this  research  effort  by  obtaining 
corresponding  probe  measurements  in  the  test  section  after  conversion  to 
solid  walls  by  the  application  of  tape  to  the  top  and  bottom  walls. 

While  reviewing  the  empty-tunnel  runs  performed  as  part  of  the 
experiments  in  Ref.  9,  however,  it  was  observed  that  although  the  uniformity 
of  the  flow  was  comparable  to  that  measured  in  earlier  experiments,  plenum 
pressure  control  valve  settings  to  achieve  uniform  flow  were  quite  different. 
In  particular,  at  Mw  ■  0.8  and  0.9,  the  first  three  upper  plenum  chambers 
were  providing  near-maximum  suction,  while  the  first  two  lower  plenum 
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chambers  were  providing  slight  blowing.  Therefore,  it  is  suspected  that 
the  flow  was  inclined  upwards  at  a  constant  angle  along  the  length  of  the 
test  section.  This  would  cause  probe  calibration  errors  and  could  lead  to 
observation  of  an  apparent  negative  cross-flow  bias  in  the  experiments. 
Unfortunately,  a  uniform  flow  inclination  along  the  length  of  the  test 
section  would  not  be  apparent  in  tho  measurements  of  differential  pressures 
across  the  static  pipe  since  those  measurements  detect  dv/dx,  rather 
than  v  itself. 

Accordingly,  in  the  most  recent  empty-tunnel  runs  at  »  0.8, 

0.9  and  0.94,  extreme  care  has  been  taken  to  insure  that  the  upper  and  lower 
valves  were  set  to  provide  comparable  amounts  of  suction  at  the  upper  and 
lower  walls.  The  asymmetry  of  the  test  section,  namely  ten  upper  plenum 
chambers  and  eight  lower  onos,  complicates  the  adjustment  of  comparable 
amounts  of  suction.  However,  it  has  been  observed  (see  Fig.  6  of  Ref.  9) 
that  relatively  large  variations  in  pipe  static  pressure  measurements 
occur  near  the  upstream  and  downstream  ends  of  the  test  section.  Thus,  by 
relaxing  tho  requirements  at  the  upstream  end  during  the  early  stages  of 
the  adjustment  process,  uniform  pressure  could  be  set  over  most  of  the  test 
section  length  with  balanced  amounts  of  suction.  Further  adjustments  then 
could  be  made  upstream  and  downstream  to  achieve  uniform  flow.  The  final 
valve  settings  do  indicate  that  comparable  suction  is  being  applied  at  the 
upper  and  lower  walls.  The  data  obtained  in  this  fashion  provided  revised 
probe  calibrations  for  zero  flow  inclination  and  were  incorporated  into  the 
data  reduction  procedure. 

q 

In  the  AEDC-sponsored  experiments  ,  the  differential  pressure 
measurements  in  uniform,  parallel  flow  indicated  an  apparent  gradient  in 
normal  velocity  at  the  upper  pipe  for  all  Mach  numbers  tested.  It  was 
concluded  that  this  gradient  is  an  artifact  of  the  pipe  construction  and 
installation.  Accordingly,  it  was  eliminated  by  treating  the  empty-tunnel 
gradient  as  a  tare  reading.  In  the  most  recent  experiments,  a  slight 
revision  of  the  empty-tunnel  gradient  tare  was  made  for  the  upper  pipe, 
and  a  small  tare  for  the  lower  pipe  was  introduced  as  well. 
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Section  4 


EXPERIMENTS  WITH  SUPERCRITICAL  FLOW  AT  THE  WALLS 
4.1  Background 

The  major  part  of  this  investigation  has  been  directed  at  flows 
which  are  supercritical  at  the  control  surfaces  and  walls.  Experiments 
were  performed  which  had  the  objective  of  demonstrating  the  segmented- 
plenum,  perforated-wall  implementation  of  the  adaptive-wall  tunnel  concept 
for  these  practically  important  flow  conditions.  In  the  previous  section, 
changes  in  the  test  section,  auxiliary  suction  system,  airfoil  model  and 
the  instrumentation  have  been  described.  All  of  these  were  accomplished 
in  order  to  provide  the  capability  for  completing  these  experiments 
successfully.  Before  summarizing  briefly  the  early  supercritical-wall 
experiments  “  and  describing  in  detail  the  most  recent  experiments  at 
"  0>9'  the  rationale  for  the  choice  of  supercritical-wall  cases  will 
be  described. 

The  lift  data  for  the  6-inch  chord  NACA  0012  airfoil,  as  tested6 
in  the  Calspan  Eight-Foot.  Tunnel  at  a  chord  Reynolds  number,  Re  ,  of 
1.00  x  10°,  showed  unusual  characteristics  at  MI)Cl  =  0.85  and  0.9,  see 
Fig.  2.  In  particular,  the  lift-curve  slope  is  very  small  at  small  cL 
and  then  has  a  gradual  break  between  =  2°  and  3°  to  a  larger  slope. 

This  behavior  was  not  observed  at  =  0.8  and  below  (see  Fig.  5  of  Ref.  6), 
or  at  »  0.925  or  0.95,  as  shown  in  Fig.  2.  The  airfoil  surface  pressure 
data  for  these  cases6'8  show  that  both  the  upper  and  lower  airfoil  surfaces 
are  supercritical  over  a  fraction  of  the  chord  length.  It  might  bo  expected 
that  the  shock-wave  formations  on  the  upper  and  lower  surfaces  would  be 
sensitive  to  wall-interference  effects,  so  these  test  conditions  promised 
to  be  realistic  tests  of  the  Calspan  adaptive -wall  implementation.  Accordingly, 
operation  at  0.85  and  0.9  was  originally  selected  for  the  supercritical - 
wall  experiments. 
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Later,  while  conducting  supercritical -wall  tests  at  M^-  0.9,  diffi¬ 
culties  encountered  in  the  iteration  process  prompted  the  reconsideration  of 
the  Eight-Foot  Tunnel  data.  A  close  examination  of  the  Eight-Foot  Tunnel  data 
at  ■  0.8S  and  0.9  for  low  oC  indicates  considerable  scatter.  This  appar¬ 
ently  results  from  relatively  large  changes  in  the  shock-wave  locations  on  the 
airfoil  due  to  small  changes  in  OC  .  Furthermore,  the  data  in  Fig.  2  and  in 
Fig.  5  of  Ref.  6  show  that  there  was  a  small  positive  flow  angularity  in  the 
Eight  Foot  Tunnel  for  all  Mw  except  0.9,  where  it  was  -0.62°.  The  intercept 
at  ■  0.9  is  suspect,  however,  because  of  the  aforementioned  scatter  in 
the  data.  This  intercept  has  been  used  in  the  experiments  presented  here, 
but  the  uncertainties  associated  with  the  Eight-Foot  Tunnel  data  in  this  Mach 
number  range  might  have  unduly  contributed  to  the  difficulties  encountered  in 
providing  definitive  verification  of  interference-free  flow.  In  principle, 
these  uncertainties  could  be  removed  by  operating  at  a  higher  Mach  number. 

In  these  studies,  however,  this  was  prevented  by  the  limitations  in  the  One- 
Foot  Tunnel  capabilities . 

One  final  aspect  of  the  data  concerns  the  Reynolds  number.  In 

the  experiments  with  the  4%-blockage  model  in  the  One-Foot  Tunnel,  Re  is 

6  ^ 
limited  to  0.67  x  10  .  This  limitation  arises  because  the  tunnel  cannot 

be  run  continuously  at  conditions  which  provide  for  Rec  ■  1.00  x  10^  with 

the  smaller  model.  Both  these  values  are  in  a  range  where  significant 

differences  in  the  detailed  and  integrated  flow  characteristics  can  be 

expected  as  a  function  of  Re.  Therefore,  a  question  remains  about  the 

pressure  distributions,  including  shock-wave  locations,  that  are  to  be 

expected  in  interference- free  flow  about  the  4%-blockage  model. 

4 . 2  Initial  Supercritical  Wall  Experiments 

The  first  supercritical-wall  experiments  with  the  4%-blockage 

O 

NACA  0012  model  were  carried  out  at  M  «  0.85  and  o6*  1°.  This  case 
is  in  the  range  where  the  lift-curve  slope  is  very  small,  see  Fig.  2. 

These  experimen. s  were  inconclusive  because  of  flow-field  unsteadiness. 

Wall  control  was  used  to  obtain  a  first  iterative  step  toward  unconfined 
flow,  but  the  shock  wave  on  the  lower  surface  fluctuated  over  about  15% 
of  the  chord.  Subsequent  attempts  to  iterate  at  this  test  condition  did  net 
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lead  to  a  steady  flow  field,  and  it  was  concluded  that  this  test  condition 
was  not  suitable  for  iteration  at  this  stage  of  the  tunnel  development. 

Next,  experiments  were  performed  with  the  4%-blockage  model  at 
M  ^  -0.8  and  CL*  4®,  for  which  there  is  no  shock  wave  on  the  lower  airfoil 
surface  and  the  upper  shock  extends  almost  to  the  static  pipe.  It  was  felt 
that  this  case  should  be  a  less  severe  test  of  the  tunnel  and  could  provide 
a  basis  for  proceeding  to  higher  M  ^  .  The  iteration  in  this  case  exhibited 
steady  flow  at  each  step  and  convergence  was  approached  as  shown  by  compari¬ 
sons  between  the  measured  data  and  evaluations  of  the  functional  relationships 
after  three  iterative  steps.  After  the  third  step,  the  relationships  were 
in  good  agreement  with  the  worst  discrepancies  occurring  in  u  at  the  lower 
control  surface  downstream  of  the  model.  A  comparison  of  the  airfoil 
pressure  distributions  also  suggested  that  the  flow  field  was  approaching 
unconfined  flow.  A  fourth  iterative  step  was  attempted  at  this  test  condi¬ 
tion.  The  u  components  were  set  close  to  the  desired  values,  but  the 
resulting  v  components  were  in  much  worse  agreement  with  the  corresponding 
unconfined- flow  distributions.  This  iterative  step  was  attempted  twice 
with  similar  results,  which  seemed  to  indicate  a  tunnel  cross-flow  condition. 
The  airfoil  pressure  distribution  confirmed  that  this  was  apparently  a 
divergent  step.  The  reasons  for  this  behavior  were  not  clear.  It  was  clear, 
however,  that  an  improved  measurement  technique  for  the  v  component  was 
necessary . 


Accordingly,  the  new  static-pipe  technique  for  determining  the  v 

9 

component  was  developed  and  this  case  was  investigated  further  .  Two 
converging  steps  in  an  experimental  iteration  were  accomplished  at  M  —  ■  0.8 
and  oC  ■  4®.  From  runs  made  prior  to  beginning  the  actual  iteration,  it 
was  concluded  that  undesired  cross-flow  conditions  upstream  of  the  model 
could  be  minimized  in  the  first  iterative  step  by  initiating  the  wall- 
control  adjustment  with  the  upstream  valves  set  at  their  tunnel -empty 
positions  at  that  M  .  These  two  iterative  steps  yielded  sufficient 
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information  to  proceed  to  a  flow  with  both  walls  supercritical.  Accordingly, 
the  same  model  was  tested  at  the  same  oL  ,  but  with  M  increased  to  0.9. 

Initial  experiments  at  ■  0.9  and  0C  »  4°  were  performed^ 
with  the  upstream  control  valves  at  their  tunnel -empty  settings  and  the 
remaining  valves  at  their  M  w  »  0.8  and  <3C=>  4°  settings.  The  measured 
components  at  the  control  surfaces  exhibited  a  reasonable  distribution 
from  the  beginning  of  the  test  section  to  the  vicinity  of  the  airfoil  trail¬ 
ing  edge  with  these  settings.  Downstream  of  this  point,  however,  the  flow 
was  choked  and  supersonic  flow  persisted  to  the  end  of  the  test  section. 

The  choking  was  relieved  by  increasing  the  suction  at  plenum  chambers  6  and  16 
just  upstream  of  the  model  (see  Fig.  1).  Further  control  adjustments  led  to 

reasonable  u  distributions  along  the  entire  length  of  the  control  surfaces. 
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Achievement  of  this  flow"  provided  the  starting  point  for  the  most  recent 
experiments  that  are  described  in  the  next  three  subsections. 

4 .3  Experiments  at  M  »  «  0.9  and  CL  »  3° 

The  angle  of  attack  of  the  4%-blockage  NACA  0012  airfoil  was 
readjusted  to  the  nominal  condition  of  0C  =  3°,  which  corresponds  to  the 
geometric  angle  of  3.62°  in  the  One-Foot  Tunnel  when  the  zero-lift  offset 
in  the  Eight-Foot  Tunnel  (see  Fig.  2)  is  taken  into  account.  Tais  case 
is  in  the  region  where  the  lift-curve  slope  is  larger  than  it  is  at  very 
small  and  provided  a  steady  case  for  which  o4  could  be  reduced  in  later 
experiments.  Adjustment  of  the  flow  at  this  condition  was  begun  with  the 
valves  in  plentun  chambers  2,  3,  14  and  15  (see  Fig.  1)  at  the  new 
empty-tunnel  settings  obtained  at  ■  0.9,  and  with  the  remaining  valves 
at  their  settings  for  the  successful  run  at  =0.9  and  0L  =  4°. 

A  first  iterative  step  toward  unconfined  flow  at  this  condition 

was  performed  by  setting  the  pressures  on  the  sides  of  the  static  pipes 

nearest  the  model  approximately  to  the  values  predicted  by  the  TSFOIL 
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computer  code.  Care  was  taken  to  insure  that  comparable  amounts  of 


suction  were  used  at  the  upper  and  lower  walls  upstream  of  the  model. 

At  the  upper  wall,  blowing  was  required  in  plenum  chambers  8,  9,  10  and  11 
in  this  and  all  subsequent  iterative  stops  at  dC  »  38,  2°  and  1°.  The 
normal  velocity  components,  as  measured  by  the  new  static-pipe  technique, 
are  presented  in  Figs.  3  and  4  at  the  upper  and  lower  control  surfaces, 
respectively.  These  measured  v  distributions  were  used  as  the  boundary 
conditions  at  the  control  surfaces  to  evaluate  the  functional  relationships 
which  must  be  satisfied  in  unconfined  flow  according  to  solutions  of  the 
transonic  small -disturbance  equations.  The  measured  streamwise  disturbance 
velocity  components,  also  found  by  using  the  new  static-pipe  technique, 
are  presented  in  Figs.  5  and  6  at  the  upper  and  lower  control  surfaces, 
respectively,  along  with  the  results  of  the  functional-relationship  evalua¬ 
tions  u[v]  .  Reasonable  agreement  is  observed  upstream  of  the  shock, 
but.  the  agreement  is  unsatisfactory  downstream  of  the  shock.  In  addition, 
the  pressure  distribution  on  the  airfoil  was  in  fair  agreement  with  the 
Eight-Foot  Tunnel  data  on  both  surfaces  up  to  the  shock  wave.  On  the 
lower  surface  of  the  airfoil,  the  shock  is  at  the  trailing  edge  as  it  was 
in  the  Eight-Foot  Tunnel  tests.  On  the  upper  surface,  however,  the  shock 
is  too  far  forward  of  the  Eight-Foot  Tunnel  position  by  about  25%  of  the 
chord  (x/c  ■  0.45  instead  of  x/c  -  0.70)  and  agreement  with  the  Eight-Foot 
Tunnel  data  is  very  poor  downstream  of  the  shock. 

A  second  iterative  step  was  then  taken  based  on  the  results  of 
the  first  step.  Small  changes  in  the  static  pressure  settings  from  the 
previous  step  were  taken,  corresponding  to  an  iterative  relaxation  factor 
of  0.1.  This  small  value  of  the  relaxation  factor  had  been  found  to  be 
necessary  in  earlier  experiments.  The  resulting  functional-relationship 
evaluations,  however,  indicated  that  very  little  improvement  had  actually 
been  accomplished.  Accordingly,  the  relaxation  factor  was  increased 
to  0.5  for  the  third  iterative  step. 
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Figure  3,  MEASURED  NORMAL  DISTURBANCE  VELOCITY  DISTRIBUTION 
-  0.9,  X  ■  3® ,  4%-BLOCKAGB  NACA  0012  AIRFOIL 
FIRST  ITERATIVE  STEP  AT  UPPER  CONTROL  SURFACE , 
h/c  -  1.0. 
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Figure  4.  MEASURED  NORMAL  DISTURBANCE  VELOCITY  DISTRIBUTION, 
M l)a  -  0.9,  .>:«  3*.  4% -BLOCKAGE  NACA  0012  AIRFOIL, 
FIRST  ITERATIVE  STEP  AT  LOWER  CONTROL  SURFACE, 
h/c  *  -1.0. 
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Figure  S.  STREAMWISB  DISTURBANCE  VBLOCITY  DISTRIBUTIONS, 

M*  -  0.9,  *>C  ■  3°,  4%-BL0CKAGE  NACA  0012  AIRFOIL, 
PIRST  ITERATIVE  STBP  AT  UPPER  CONTROL  SURFACE, 
h/c  "1.0 
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The  agreement  of  the  measured  u  distributions  at  the  control 
surfaces  in  the  third  iterative  step  with  the  distributions  u  [v]  from 
the  functional -relationship  evaluations  improved  significantly  downstream 
of  the  model  and  remained  good  elsewhere  as  shown  in  Pigs.  7  and  8.  The 
shock  wave  moved  rearward  on  the  upper  surface  of  the  airfoil  and  is  now 
too  far  forward  by  about  15%  of  the  chord  (x/c  -  0.55  instead  of  x/c  «  0.70). 

Overall,  the  third  iterative  step  was  a  significant  improvement  and 
showed  that  sufficient  control  was  available  to  iterate  toward  unconfined 
flow.  Moreover,  the  new  static-pipe  technique  is  a  great  improvement  over 
the  use  of  the  flow-angle  probes  alone,  especially  with  respect  to  ease 
and  repeatability  of  the  measurements.  Nevertheless,  certain  limitations 
in  the  present  implementation  of  the  new  technique  were  found  to  exist  and 
certain  fundamental  questions  about  the  pipe  characteristics  remain  to  be 
resolved.  These  will  become  clearer  by  discussing  the  representative  results 
that  are  shown  in  Figs.  3  and  4. 

The  most  obvious  characteristic  in  the  results  of  Figs.  3  and  4 
is  the  dashed  line  representing  the  break  in  the  curves  in  the  vicinity  of 
the  shock  wave.  There  are  three  problems  near  and  downstream  of  the  shock. 
First,  the  pipe  orifices  are  spaced  1-inch  apart,  so  that  there  is  a  gap 

g 

of  A  (x/c)  «0.25  between  measurements.  Second,  the  slender-body  theory 
relating  pipe  pressures  to  u  and  dv/dx  is  not  strictly  applicable  when 
the  disturbance  causes  formation  of  a  supersonic  pocket,  followed  by  a 
shock  such  as  occurs  at  both  pipes  in  this  case.  Third,  the  response  of 
the  pipe  boundary  layer  to  the  shock  wave  is  not  known.  Proper  assessment 
of  the  inviscid  shock-pipe  interaction  and  the  possible  effects  on  the  pipe 
boundary  layer  were  beyond  the  scope  of  this  investigation,  as  was  the 
development  of  a  pipe  traversing  mechanism.  In  order  to  overcome  these 
problems,  it  has  been  assumed  that  the  slender-body  analysis  applies 
independently  upstream  and  downstream  of  the  shock.  Then  the  curves  on 
each  side  of  the  shock  are  extrupolated  graphically  as  shown  by  the 
dashed  lines. 
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Figure  7.  STREAMWISE  DISTURBANCE  VELOCITY  DISTRIBUTIONS, 

■  0.9,  oC  .  3*.  4%-BLOCKAGE  NACA  0012  AIRFOIL 
THIRD  ITERATIVE  STEP  AT  UPPER  CONTROL  SURFACE , 
h/c  *  1.0. 
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Figure  8.  STREAMWISE  DISTURBANCE  VBLOCITY  DISTRIBUTIONS, 

■  0.9,  0t«  3°,  4%- BLOCKAGE  NACA  0012  AIRFOIL 
THIRD  ITERATIVE  STEP  AT  LOWER  CONTROL  SURFACE, 


In  each  of  Pigs.  3  and  4,  the  data  for  dv/dx  up  to  the  shock  were 
fit  by  cubic  splines  and  integrated  up  to  the  last  data  point  before  the 
shock  with  the  constant  of  integration  found  by  a  least-squares  fit  to 
the  probe  data  over  that  interval.  Figures  3  and  4  are  typical  of  all  the 
results  upstream  of  the  shock  in  that  the  generally  smooth  v  distributions 
are  similar  to  those  expected  from  theoretical  considerations.  Downstream 
of  the  shock,  a  similar  spline  fit  and  integration  were  performed.  At  the 
lower  pipe,  the  probe  measurement  at  x/c  =  0.8275  was  not  considered  in  the 
fit  to  determine  the  constant  of  integration  because  its  location  with 
respect  to  the  shock  was  not  known  precisely,  nor  was  its  accuracy. 

The  results  in  Figs.  3  and  4  downstream  of  the  shock  are  typical 
of  most  of  the  experiments  to  be  described  below,  especially  with  regard 
to  the  variations  in  v  downstream  of  the  shocks .  These  unusual  variations 
give  rise,  of  course,  to  the  variations  in  u[v]  shown  in  Figs.  5  to  8. 

In  particular,  the  asymptotic  approach  to  v  -  0  far  downstream  was  not 

always  from  negative  values  of  v  at  the  upper  pipe  as  it  is  in  Fig.  3, 
while  in  most  cases,  the  approach  at  the  lower  pipe  was  from  the  negative 
side  as  shown  in  Fig.  4.  Global  considerations  for  the  disturbances  gen¬ 
erated  by  the  model  in  unconfined  flow  indicate  that  v  ■  0  should  be 
approached  from  the  negative  side  for  the  upper  pipe  and  from  the  positive 
side  for  the  lower  pipe.  There  are  several  possible  reasons  for  deviations 
from  this  expected  behavior.  First,  the  applicability  of  the  subsonic, 
shock-free  theory  relating  the  pipe  pressures  to  u  and  dv/dx  is  in  doubt 
near  the  shock.  Second,  the  flow  control  during  early  steps  in  the  adaptive 
wall  iterative  scheme  may  be  such  that  the  flow  is  so  far  from  interference - 
free  conditions  that  the  measurements  are  actually  correct  despite  their 
unexpected  behavior.  Finally,  the  probe  measurements  used  to  fix  the 
constant  of  integration  may  be  in  error,  especially  at  the  lower  pipe, 
where  only  the  single  measurement,  at  x/c  ■  2.515  is  available. 
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These  lingering  uncertainties  about  the  interpretation  of  the 
static-pipe  measurements  led  to  the  conclusion  that  the  third  iterative  step 
was  probably  at  the  limits  of  measurement  accuracy,  so  far  as  further 
iterations  were  concerned.  Therefore,  it  was  decided  instead  to  reduce  oi 
successively  to  2®  and  Is,  so  that  an  investigation  could  be  made  of  these 
cases  where  the  lift-curve  slope  is  smaller.  It  was  in  this  range  at 
M  p  m  0.85  that  steady  flow  could  not  be  achieved. 

Before  these  lower  CL  cases  were  investigated,  however,  several 
runs  were  made  at  M  w  «0.9  and  (X.  ■  3°  to  examine  the  effects  of 
increased  Reynolds  number.  In  order  to  achieve  these  conditions,  the 
tunnel  pressure  conditions  were  changed.  However,  no  valve  setting  changes 
were  made  in  the  auxiliary  suction  and  pressure  systems.  The  power  require¬ 
ments  and  heating  behavior  of  the  tunnel  were  such  that  run  time  was  limited 

to  about  20  minutes  at  the  higher  Re  .  This  time  did  permit  acquisition  of 

c 

pressure  data  on  the  upper  surface  of  the  airfoil. 

As  Rec  was  increased  from  0,67  x  106  to  0.83  x  106,  on  the 

upper  airfoil  surface  decreased  by  about  30%  (acceleration  of  the  flow) 

everywhere  upstream  of  the  shock  wave,  which  moved  rearward  by  about  2.5% 

of  the  chord.  As  Re  was  increased  further  to  1.00  x  106,  C  decreased 

c  p 

by  an  additional  30%  ahead  of  the  shock,  which  moved  forward  by  about  7.5% 
of  the  chord  so  that  it  was  5%  ahead  of  its  location  at  Rec  «  0.67  x  106. 

The  shock  became  stronger  as  Re^  increased.  This  behavior  apparently  results 
from  the  interaction  of  the  revised  pressure  levels  in  the  tunnel  with  those 
in  the  auxiliary  suction  and  pressure  systems.  That  is,  wall-interference 
effects  are  introduced  by  the  mismatch  in  pressures.  Complete  readjustment 
of  the  wall  control  at  each  flow  condition  would  be  necessary  to  assess  the 
effect  of  Rec  on  the  model  properly.  This  was  not  feasible  with  the  limited 
run  time.  These  tests  definitely  established  the  sensitivity  of  the  flow 
to  Reynolds  number  and  indicated  the  necessity  to  repeat  the  iteration 
procedure  whenever  M  ^  ,  tX  or  Rec  are  changed. 
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The  airfoil  angle  of  attack  was  reduced  to  the  nominal  value 
(Fig.  2)  of  (C  “  2°,  which  corresponds  to  the  geometric  angle  of  oC  ■  2.62° 
in  the  One-Foot  Tunnel.  The  first  step  in  the  iterative  procedure  was  made 
by  retaining  the  plenum  pressure  control  valves  at  their  final  settings  for 
the  OC  m  3*  iterations  and  simply  acquiring  a  full  set  of  static- pipe  and 
probe  data.  The  agreement  at  both  control  surfaces  between  the  measured  u 
distributions  and  the  u[v]  obtained  from  the  functional -relationship 
evaluations  is  generally  comparable  to  that  shown  in  Figs.  S  and  6  for  the 
first  iterative  step  at  OC  =  3°.  The  largest  differences  still  are  down¬ 
stream  of  the  shock  waves .  The  shock  on  the  upper  surface  of  the  airfoil 
is  at  x/c  ■  0,55,  as  compared  to  x/c  =  0.75  in  the  Eight-Foot  Tunnel  data. 

On  the  lower  surface,  the  shock  is  at  x/c  *  0.95  instead  of  0.975  as  in  the 
Eight-Foot  Tunnel  data. 

A  second  iterative  step  was  made  with  a  relaxation  factor  of  0.5, 
as  had  proved  better  at  OS-  -  3°.  This  choice  was  satisfactory  here,  too, 
inasmuch  as  the  agreement  at  this  step  is  comparable  to  that  shown  in 
Figs.  7  and  S  for  the  third  step  at  <*  3°.  The  shock  on  the  upper  surface 
of  the  airfoil  moved  rearward  to  x/c  ■  0.65,  while  the  lower  surface  shock 
moved  slightly  rearward. 

Overall,  then,  these  results  at  CL-  2°  reflected  the  behavior  at 
OL  -  3#.  Hie  flow  was  again  steady,  without  evidence  of  shock-wave  fluctu¬ 
ations.  At  the  lower  control  surface,  the  measured  v  distributions  at 
both  iterative  steps  were  very  similar  to  those  at  OC  «-•  3®,  (see  Fig.  4). 

The  asymptotic  approach  to  zero  downstream  was  still  contrary  to  that 
expected  for  interference- free  flow.  At  the  upper  control  surface  in  both 
iterative  steps,  however,  there  was  a  departure  from  the  behavior  shown  in 
Fig.  3;  that  is,  the  asymptotic  approach  to  zero  downstream  was  from  the 
positive  side.  Despite  these  questions,  the  iterations  for  this  case  were 
terminated  at  the  second  step  so  that  the  OC  ■  1°  case  could  be  investigated 
more  thoroughly. 


4.5  Experiments  at  Mm»  g  0.9  and  CL  •  1° 

The  angle  of  attack  was  reduced  to  the  nominal  setting  of  0LU  1°, 
which  corresponds  to  the  geometrical  angle  of  1.62°  in  the  One-Foot  Tunnel. 
The  first  and  second  iterative  steps  were  carried  out  as  they  had  been  for 
CL  a  2°.  That  is,  the  first  step  was  made  with  the  plenum  control  valves 
at  their  final  settings  for  CL  *  2®  and  the  second  step  was  taken  with  an 
iterative  relaxation  factor  of  0.5.  The  basic  trend  of  the  results  for 
the  distributions  is  similar  to  the  first  and  second  steps  at  DC  *  2° 
and  the  first  and  third  steps  at  0L  =  3° .  However,  the  shock  wave  on  the 
upper  surface  of  the  airfoil  moved  from  x/c  =  0.65  at  the  first  step  to 
x/c  ■  0.85  at  the  second.  This  is  aft  of  its  location  at  x/c  ■  0.725  in 
the  Eight-Foot  Tunnel  tests.  At  the  lower  surface  of  the  airfoil,  the 
shock  was  slightly  forward  of  its  Eight-Foot  Tunnel  location  at  the  first 
step,  but  subsequently  moved  aft.  The  flow  was  steady  at  each  step  and 
there  was  no  exceptional  difficulty  in  setting  the  second  step. 

A  third  iterative  step  was  taken,  again  with  a  relaxation  factor 
of  0.5.  The  resulting  measu^id  v  components  at  the  control  surfaces  are 
given  in  Figs.  9  and  10.  These  are  used,  in  turn,  as  the  boundary  conditions 
for  the  functional -relationship  evaluations  u[v],  results  of  which  are 
compared  with  the  measured  u  data  in  Figs.  11  and  12.  The  agreement  in  the 
u  distributions  at  both  control  surfaces  is  far  superior  to  that  observed 
for  any  of  the  previous  cases.  However,  there  are  still  some  important 
differences  downstream  of  the  shocks.  The  chordwise  distributions  of  the 
airfoil  surface  pressures  for  the  third  iterative  step,  presented  in  terms 
of  ftCp  ,  are  compared  with  the  corresponding  Eight-Foot  Tunnel  data  in 
Fig.  13.  The  shock  waves  on  both  surfaces  are  nearly  coincident  for  the 
two  flows.  However,  there  is  an  appreciable  deviation  of  the  pressures  on 
the  upper  surface  ahead  of  the  shock  and  the  One-Foot  Tunnel  data  appear 
to  have  a  smoother  variation  than  the  Eight-Foot  Tunnel  data.  The  reason 
for  these  differences  is  not  known. 
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Figure  9.  MEASURED  NORMAL  DISTURBANCE  VELOCITY  DISTRIBUTION, 
M*,  -  0.9,  OC-  3°,  4% -BLOCKAGE  NACA  0012  AIRFOIL, 
THIRD  ITERATIVE  STEP  AT  UPPER  CONTROL  SURFACE, 
h/c  ■  1.0. 
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Figure  10.  MEASURED  NORMAL  DISTURBANCE  VBLOCITY  DISTRIBUTION, 
M^  *  0.9,  -  1°,  4%-BLOCKAGE  NACA  0012  AIRFOIL 

THIRD  ITERATIVE  STEP  AT  LOWER  CONTROL  SURFACE. 
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Figure  12.  STRBAMWISB  DISTURBANCE  VELOCITY  DISTRIBUTIONS, 

Mo*  «  0.9,  OL  m  1*,  4 BLOCKAGE  NACA  0012  AIRFOIL 
THIRD  ITERATIVE  STBP  AT  LOWER  CONTROL  SURFACE, 
h/c  -  -1.0. 


37 


The  v  distributions  measured  at  the  third  iterative  step  follow 
the  trend  expected  from  theoretical  predictions  for  unconfined  flow.  That 
is,  the  nonuniform  variation  in  the  distribution  downstream  of  the  shock 
at  the  upper  control  surface  (Fig.  9)  is  not  as  severe  as  it  was  for  the 
first  iterative  step  at  06  »  3"  (Fig.  3)  and  for  most  of  the  intervening 
cases.  Also,  the  asymptotic  approach  to  v  *  0  is  satisfactory.  At  the 
lower  control  surface  (Fig.  10),  too,  the  v  distribution  is  more  reasonable 
with  only  a  small  overshoot  in  the  asymptotic  behavior  downstream.  This 
overshoot  may  result  from  uncertainty  in  the  probe  measurement  at  x/c  ■  2.515. 

Overall,  the  trend  of  the  iterations  for  all  three  angles  of 
attack  at  ■  0.9  is  towards  convergence.  Improvements  were  made  in  the 
agreement  between  the  measured  u  distributions  and  the  u[v]  computed  by 
the  functional -relationship  evaluations  at  each  iterative  step.  Thus,  it 
would  appear  that  the  entire  angle-of-attack  range  can  be  iterated  to  a 
reasonable  approximation  of  interference- free  flow  conditions.  Further 
refinement  of  the  results  achieved  here  should  be  possible  after  resolution 
of  the  outstanding  questions  about  the  static-pipe  data  and  the  Reynolds 
number  effects, 
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Section  5 


CONCLUDING  REMARKS 

Q  i 

The  new  static-pipe  technique  for  determining  the  streamwise  and 
normal  disturbance  velocity  components  provided  significant  improvements 
in  the  ease  and  reliability  of  the  normal  component  measurements.  The 
extensive  use  of  the  technique  that  has  been  accomplished  herein,  however, 
has  led  to  uncertainties  about  the  interpretation  of  the  pipe  measurements 
in  some  regions  of  the  flow  field.  In  particular,  results  in  the  vicinity 
of  the  interaction  of  a  shock  wave  from  the  model  with  the  pipe  has  raised 
several  questions  which  must  be  answered  before  additional  use  of  the 
technique  is  undertaken.  Although  the  new  static  pipes  have  greatly  improved 
the  definition  of  the  normal  velocity  distribution  near  the  model,  an  even 
more  refined  definition  is  needed  in  regions  of  rapid  gradients,  especially 
near  a  shock,  A  translation  mechanism  for  the  static  pipes  could  provide 
this  capability.  The  subsonic  slender-body  theory  developed  in  Ref.  9 
should  be  replaced  by  a  transonic  theory,  particularly  in  the  vicinity  of 
the  shock,  in  order  to  provide  an  accurate  relationship  between  the  pressure 
measurements  on  the  pipe  and  the  two  velocity  components.  Furthermore,  an 
analysis  of  the  effects  of  the  pipe  boundary- layer  on  4 he  relationship 
between  the  pressures  and  velocity  components  should  be  performed.  Finally, 
experiments  should  be  made  in  which  the  results  of  the  static-pipe  measure¬ 
ments  can  be  compared  with  independent  measurements  of  the  two  velocity 
components.  Successful  resolution  of  these  questions  could  lead  to  use  of 
the  new  static-pipe  technique  for  measurements  in  flow  fields  other  than 
those  in  adaptive-wall  wind  tunnels. 

Successful  experimental  iterations  toward  interference-free  flow 
have  been  achieved  at  a  free-stream  Mach  number  of  0.9  and  angles  of  attack 
of  3*,  2*  and  Is.  As  the  iterations  proceeded,  significant  improvements 
were  observed  in  the  agreement  between  the  measured  streamwise  disturbance 
velocity  distributions  and  the  corresponding  distributions  resulting  from 
the  exterior- flow  functional-relationship  evaluations.  These  evaluations 
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used  the  measured  normal  velocity  distributions  as  boundary  conditions 
for  transonic  small -disturbance  calculations.  Therefore,  we  conclude  that 
the  existing  perforated-wall,  segmented-plenum  implementation  of  the 
adaptive-wall  concept  can  provide  sufficient  control  to  iterate  when  locally- 
supersonic  flow  and  shock  waves  extend  to  both  tunnel  walls.  In  the  wall- 
control  adjustment  procedure,  care  must  be  exercised  to  balance  the  amount 
of  suction  in  the  upper  and  lower  plenum  chambers  upstream  of  the  model. 

This  is  necessary  to  avoid  introducing  unwanted  flow  inclination  upstream 
of  the  model.  Furthermore,  positive  pressure  was  necessary  in  the  upper 
plenum  chambers  downstream  of  the  model  in  order  to  provide  the  blowing 
velocity  through  the  wall  from  the  plenum  chambers  into  the  tunnel  that 
was  required  to  achieve  control  of  the  flow. 

We  believe  that  an  improved  two-dimensional  test  section  with  the 
same  perforated-wall,  segmented-plenum,  wall-control  technique  is  possible. 
Those  improvements,  which  are  described  in  Appendix  II,  should  provide  more 
refined  control  in  the  near  vicinity  of  the  model  for  supercritical-wall 
cases.  The  refinements  would  be  accomplished  by  changing  plenum- chamber 
size,  particularly  by  reducing  the  length  of  those  chambers  near  the  model 
where  the  gradients  in  the  disturbance  velocity  components  are  greatest. 

Also,  walls  with  different  open-area  ratios  for  smaller  losses  in  the 
auxiliary  suction  system  might  be  possible  away  from  the  immediate  vicinity 
of  the  shock  waves.  Finally,  symmetrical  plenum  segmentation  on  the  upper 
and  lower  walls  would  provide  greater  flexibility  of  operation,  particularly 
for  calibration  and  routine  checks  of  the  instrumentation,  both  without  a 
model  present  and  with  a  symmetrical  model  at  zero  angle  of  attack. 

The  numerical  simulation  investigation  is  presented  in  detail  in 
the  AIAA  paper  that  is  attached  to  this  report,  and  extensive  concluding 
remarks  are  given  there.  However,  it  is  worthwhile  to  summarize  them  here. 
The  low-speed  numerical -simulation  methodology  gives  a  realistic  description 
of  the  behavior  of  a  perforated- wal 1 ,  segmented-plenum  test  section,  as 


demonstrated  by  baseline  numerical  simulations  of  a  representative  case. 
Influence  functions  calculated  for  changes  in  the  velocity  components  at 
the  control  surfaces  due  to  a  change  in  the  pressure  in  one  plenum  chamber 
at  a  time,  while  holding  the  other  chamber  pressures  constant,  can  be 
approximated  reasonably  well  by  universal  curves  based  on  linearized  theory. 
Extension  of  the  numerical  simulation  technology  to  flows  which  are  super¬ 
critical  at  the  walls  could  be  used  as  a  basis  for  wide-ranging  studies  of 
test-section  design  and  operational  procedures.  In  particular,  simulations 
would  enable  the  development  of  procedures  for  automated  wall  adjustment, 
which  is  crucial  to  the  ultimate  practical  application  of  the  adaptive-wall 
concept.  This  is  clearly  apparent  in  light  of  the  great  care  that  was 
necessary  in  the  manual  adjustments  required  during  the  iteration  experiments 
described  in  this  report. 

As  a  result  of  our  entire  experience  with  the  two-dimensional 
test  section,  through  both  experiment  and  numerical  simulation,  and  the 
preliminary  design  studies  described  in  Appendix  II,  wo  believe  that  the 
perforated-wall,  segmentod-plenum  implementation  is  viable  for  a  fully 
three-dimensional  test  section.  This  would  be  especially  so  if  additional 
design  studies  were  to  show  that  adequate  control  would  be  possible  with 
even  a  smaller  number  of  plenum  chambers  than  the  64  discussed  in  Appendix  II. 
Appropriate  adaptations  of  the  numerical  simulation  methodology  could  be  used 
to  perform  these  design  studies .  The  new  static-pipe  technique  has  been 
g 

shown  theoretically  to  be  applicable  to  the  measurement  of  three- 
dimensional  disturbance  velocity  fields  when  suitably  modified.  This 
technique  is  an  attractive  possibility  for  use  in  a  three-dimensional 
adoptive-wall  test  section. 
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In  addition,  Calspan  personnel  participated  in  Adaptive  Transonic 
Wall  Meetings  held  at  ABDC  on  19  January  1978,  10  August  1978,  19  January  1979 , 
12  July  1979,  26  February  1980  and  16  January  1981.  Research  performed  under 
this  contract  was  presented  at  these  meetings. 
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THREE-DIMENSIONAL  ADAPTIVE-WALL  DESIGN  STUDIES 


The  objectives  of  these  studies  were  to  establish  requirements  for 
three-dimensional  adaptive-wall  test  sections  and  then  to  examine  the  adapta¬ 
bility  of  the  Calspan  One-Foot  Tunnel  to  meet  these  requirements.  The  test 
section  requirements  were  investigated  by  calculating  the  unconfined  flow 
about  a  few  representative  configurations.  Estimates  were  based  on  both 
Prandtl-Glauort  theory  and  finite-difference  solutions  of  the  tiansonic  small- 
disturbance  equations,  using  the  Bailey-Ballhaus  code;0  The  emphasis  in  the 
adaptation  study  was  on  working  toward  the  design  of  a  suitable  test  section, 
as  well  as  a  corresponding  model  configuration  with  which  to  carry  out  a 
demonstration  of  the  concept. 

Choice  of  a  suitable  model  configuration  for  demonstration  testing 
is  very  important.  The  model  should  not  only  be  representative  of  practical 
configurations,  but  should  also  have  a  large  amount  of  interference  in  a 
passive-wall  tunnel  so  that  a  clear-cut  demonstration  can  be  made.  A  wing- 
body-horizontal  tail  model  is  a  minimum  requirement  so  that  the  contribution 
of  the  horizontal  tail  to  the  pitching  moment  can  be  assessed.  This  contri¬ 
bution  should  be  sensitive  to  the  capability  of  the  chosen  adaptive-wall 
implementation  for  assuring  that  the  trailing  vortex  system  from  the  wing 

retains  its  unconfined-flow  characteristics  and  location.  The  most  suitable 

23 

existing  model  is  the  AEDC  Wall -Interference  Model  ,  which  has  a  slender 
axisymmetric  body  with  a  swept,  untapered,  const ant -thickness  wing  and  horizon 
tal  tail.  This  model  is  attractive  for  a  three-dimensional  demonstration 
because  the  wall  interference  should  be  relatively  large. 
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Bailey,  F.R.  and  Ballhaus,  W.F.  "Comparison  of  Computed  and  Experimental 
Pressures  for  Transonic  Flows  About  Isolated  Wings  and  Wing-Fuselage 
Configurations"  NASA-SP-34?-Part  2,  March  1975. 


Before  the  above  model  was  considered,  initial  calculations  of  the 
unconfined  flows  about  two  different  wing  configurations  were  made  with  the 
Bailey-Ballhaus  code.  One  wing  was  swept,  tapered  and  untwisted  while  the 
other  was  rectangular  in  planform  and  untwisted.  Thus,  their  flowfields, 
which  were  calculated  for  a  single  angle  of  attack,  each  at  a  Mach  number 
of  0.9,  were  quite  distinctive.  A  computer  program  was  written  to  evaluate 
the  disturbance  velocity  components  at  walls  or  control  surfaces  of  rectangu¬ 
lar  cross  section.  This  program  reads  the  velocity  potential  for  unconfined 
flow  throughout  the  field,  as  saved  on  tape  in  the  Bailey-Ballhaus  program, 
and  differentiates  and  interpolates  numerically  to  get  the  desired  flowfield 
velocities.  It  also  estimates  the  separate  effects  of  thickness  and  lift  so 
that  superposition  can  be  used  to  approximate  the  velocity  field  for  other 

values  of  these  parameters.  Superposition  is  not  strictly  valid  for  nonlinear 
supercritical  flows.  However,  it  does  provide  useful  estimates  of  the 

velocity  magnitudes  and  distributions,  especially  when  the  flow  at  the  field 
location  is  subcritical . 

Prandtl-Glauert  calculations  based  on  line  distributions  of 
thickness  and  lift  for  rectangular  planforms  were  useful  for  estimating  the 
effects  of  free-stream  Mach  number  on  the  disturbance  velocity  distributions. 
In  both  two  and  three  dimensions,  an  increase  in  Mach  number  shifts  closer 
together  the  maxima  and  minima  of  the  normal  and  streamwise  velocity  distribu¬ 
tions  on  the  upper  and  lower  walls.  For  bodies,  slender-body  theory  was  used 
to  obtain  an  axial  source  distribution  for  a  given  axiiymmetric  body-cross¬ 
section  distribution.  The  source  distribution  is  then  approximated  by 
discrete  source  elements,  which  are  used  in  conjunction  with  Prandtl-Glauert 
theory  to  calculate  the  velocity  components  at  the  control  surface  and  wall 
locations . 

Once  these  three-dimensional  unconfined  flow  fields  had  been 
obtained,  we  began  to  consider  test  section  designs  so  that  we  could  define 
more  fully  the  outstanding  problems  which  require  resolution  before  a 
definitive  design  can  be  achieved,  licrly  in  the  investigation,  we  examined 
the  prospects  for  a  reflection-plane  test  section  with  control  on  only  three 


walls.  This  is  attractive,  principally,  because  a  larger  model  can  be  accom¬ 
modated  for  a  given  tunnel  cross-sectional  area,  thus  providing  more  uncon¬ 
trolled  wall  interference  to  be  eliminated  in  the  demonstration.  However,  a 
convincing  demonstration  would  require  that  the  same  model  also  be  tested  in 
a  larger  tunnel  to  obtain  substantially  irterference-free  data  for  comparisons. 
It  would  be  very  difficult  to  account  properly  for  the  differences  in  the 
reflection-plane  boundary  layers  on  tests  in  two,  very  different  tunnel  sizes. 
Moreover,  for  model  configurations  that  are  symmetrical  laterally,  we  believe 
that  it  is  important  to  demonstrate  that  the  walls  can  be  adjusted  so  as  to 
insure  the  same  symmetry  in  the  flow  field.  Also,  demonstrations  should  be 
made  for  moderately  non-symmetrical  model  attitudes.  Therefore,  we  concluded 
that  a  fully  three-dimensional  test  section  should  be  built  with  active 
control  on  all  four  walls. 

We  decided  that  an  extension  of  our  two-dimensional  design 
principles  would  be  used  as  a  base  line.  That  is,  the  test  section  would  have 
perforated  walls,  behind  which  would  be  segmented  plenum  chambers  with  individ¬ 
ual  pressure  control.  Accordingly,  the  next  step  in  the  design  process  was 
to  reconsider  the  present  two-dimensional  design  in  light  of  our  operational 
experience.  It  appears  that  smaller  plenum  chambers  near  the  model  would  be 
advantageous  to  accommodate  the  rapid  variations  that  occur  in  the  velocity 
components  along  and  normal  to  the  walls  at  high  Mach  numbers.  The  original 
design  was  based  on  distributions  calculated  at  low  speeds.  As  mentioned 
above,  at  higher  Mach  numbers,  the  maxima  and  minima  are  closer  together  so 
that  the  distributions  of  the  velocity  components  are  no  longer  approximately 
linear  over  the  plena  near  the  model.  On  the  other  hand,  farther  away  from 
the  model,  it  appears  that  larger  plenum  chambers  might  be  possible.  In 
addition,  the  use  of  walls  with  a  larger  open-area  ratio  is  attractive  for 
reducing  the  suction  requirements  in  regions  where  no  shock  waves  are  expected. 
Where  shocks  are  expected,  the  present  22.5%  open-area  ratio  probably  should 
be  retained  because  of  its  superior  shock -cancellation  properties.  Although 
these  features  appear  attractive,  based  on  our  experience  in  two  dimensions, 
more  evidence  of  their  effectiveness  is  required  and  consequently,  further 
work  should  be  carried  out  to  resolve  these  matters. 


An  important  decision  on  the  test  section  configuration  is  the 
cross-section  shape.  Most  modem  tunnels  are  square  for  maximum  flexibility 
in  testing  a  wide  variety  of  configurations,  including  vehicles  at  very  large 
angles  of  attack.  For  the  purposes  of  an  adaptive-wall  demonstration,  we 
initially  considered  rectangular  cross  sections  with  height -to -width  ratios 
from  about  1/2  to  1 .  A  ratio  of  0.7  gives  a  good  balance  between  the  peak 
magnitudes  of  the  unconfined-flow  velocity  components  on  the  side,  top  and 
bottom  walls  for  as  large  a  model  span  as  is  practical  to  test  in  a  given 
tunnel  width.  This  still  permits  configurations  to  be  tested  at  moderate 
angles  of  attack!  This  ratio  is  also  compatible  with  the  existing  contraction 
section  of  the  Calspan  One-Foot  Tunnel. 


A  preliminary  layout  was  made  of  the  plenum  segmentation  required 
in  the  streamwise  direction,  using  the  features  set  forth  above,  based  on 
our  two-dimensional  experience.  It  was  found  that  both  the  rectangular  and 
swept  wings,  which  we  considered,  could  be  accommodated  by  the  same  basic  seg 
mentation,  provided  that  each  model  could  be  mounted  in  a  different  location 
with  respect  to  the  plenum  chambers.  We  believe  that  we  can  achieve  wall 
control  using  the  same  number  of  streamwise  plenum  chambers  as  in  two  dimen¬ 
sions,  including  provisions  for  the  fuselage-induced  disturbance  velocities. 
In  the  resulting  streamwise  breakdown,  the  velocity  components  vary  approxi¬ 
mately  linearly  over  the  extent  of  each  plenum,  A  corresponding  segmenta¬ 
tion  of  the  side,  top  and  bottom  walls  was  also  chosen  on  the  same  basis  of 
approximately  linear  variations.  Overall,  this  initial  design  had  64  indi¬ 
vidual  plenum  chambers,  with  8  laterally  at  each  of  8  streamwise  segments. 


After  completion  of  the  preliminary  layout,  alternative  tunnel 
cross-section  shapes  were  considered.  For  rectangular  cross  sections  at  a 
given  axial  station,  the  Induced  normal  velocity  for  unconfined  flow  over 
representative  wings  varies  considerably  over  the  walls  and  becomes  small  in 
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the  comers.  If  the  corners  are  modified  by  wall  segments  that  make  the 
cross  section  octagonal,  it  was  found  that  the  normal  velocity  is  nearly 
constant  over  these  additional  segments,  if  the  segments  are  proportioned 
suitably.  Also,  the  streamwise  variation  of  the  normal  velocity  is  essentially 
the  same  for  both  the  octagonal  and  rectangular  shapes .  Consideration  of 
flow  control  requirements  then  indicates  that  the  octagonal  section  does  not 
require  any  more  plenum  chambers  than  the  rectangular.  Consequently,  octagonal 
cross  sections  offer  promise  of  more  effective  flow  control,  along  with  model 
span  maximization  for  a  given  tunnel  mass  flow. 


Finally,  c  conceptual  design  for  a  new  plenum  chamber  and  header 
configuration  for  providing  both  suction  and  blowing  capability  was  developed. 
This  design  promises  a  considerable  reduction  in  the  complexity  of  the  auxi¬ 
liary  pressure  and  suction  circuits,  coitqpared  to  our  two-dimensional  design. 
However,  it  is  regarded  only  as  part  of  a  first  step  toward  a  three-dimensional 
design.  Overall,  the  conceptual  design  holds  considerable  promise  for  a  three- 
dimensional  demonstration  experiment.  Further  research  would  be  required  to 
refine  the  requirements  and  to  prepare  a  final  detailed  design. 
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Abstract 


Flow  within  the  tunnel  is  simulated  by  model¬ 
ing  the  incompressible  interaction  of  the  trans¬ 
pired  turbulent  boundary  layers  on  the  walls  with 
the  flow  over  the  airfoil.  Despite  the  fact  that 
a  finite  number  of  plenum  chambers  can  exert  only 
imperfect  control  over  the  flow,  it  Is  demonstrated 
that  one  can  still  achieve  what  la  for  all  practi¬ 
cal  purposes  unconfined  flow  about  the  airfoil. 
Velocity  differences  produced  at  control  surfacos 
outside  the  boundary  layors  by  changing  tho 
pressuro  in  ono  plenum  chnmhor  at  u  time,  holding 
the  other  chamber  prossuros  constunt,  are  presented 
as  influence  functions.  Implications  of  tiio  anal¬ 
ysis  on  tunnol  design  and  automation  tiro  described, 


NomcncJutur0 

C(,Cf  friction  coefficients  with  and  without 
*  transpiration,  respectively 

Cm  prossuro  coefficient  on  airfoil  surface, 

</*  -  >  '4- 

pressure  coefficient  In  plenum  n  , 

/*«  .  »q.  (2) 

c  airfoil  chord,  I'lg.  1 

H  boundary- layov  shapo  factor,  6/9 

h  distunco  from  model  to  tunnol  walls, 

I'ig.  1 

Ia,Iv  Influence  functions  for  a  and  v  deter¬ 
mined  from  simulations,  defined  In  Hqs. 

.  „  (5)  and  (o) 

Ju.  i  ^ir  Idotilised  Influence  functions,  defined 
In  liqs.  (ft- IP)  and  (R-20) 
in  length  of  plenum  chnmhor  n 

M  Mach  number 

p  static  prossuro 

i  dyimmic  prossuro 

U  total  velocity  component  In  31  direction 

u,v  perturbation  velocity  components  in 
X  ,  y  directions 

Au,iv  change  tn  a  ,  v  uf tor  changing 
prossuro  In  ono  plenum 
6 a  -  A  u  t-  ** ) 

vt  effective  Invlscid  norma)  velocity  at 
wall  due  to  boundary-layer  displacement 
(positive  out  of  the  wall),  lit).  (1) 
transpiration  velocity  at  wall  (positive 
out  of  tho  wall) 

X ,  coordinate  system  with  origin  at  airfoil 
quarter  chord,  Fig.  1 
stroumwlso  locution  of  center  of  plenum 
chamber  >» 

X  (  X  -  )  /  J„ 

Ue  distance  from  tunnel  contor  line  to 
control  surfaces,  Fig,  1 

i  i/h 

a  airfoil  onglo  of  attack 
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6  boundary- layer  thickness 

6  boundary -layer  displacement  thickness 

8  boundary- layer  momentum  thickness 

f,  coordinate  system  with  origin  at  the 

center  of  plenum  chamber  n  midway 
between  the  walls 

^  disturbance  velocity  potential  in 

Appendix  B 


Subscripts 

«■  ovaluated  at  control  surface  adjacent 
to  plenum  chumber  n 

e  condition  at  odge  of  boundary  layer 

n  condition  in  plenum  n 

o  evaluated  at  control  surface  opposite 

to  plonum  chamber  n 

.  oo  free- stream  condition 


1 ,  Introduction 

Tho  concept  of  an  uduptivo-wul 1  wind  tunnel 
bus  attracted  much  attention  recently1 *14  because 
of  Its  promise  for  significant  reductions  In  tunnol 
wall  interference,  particularly  in  the  transonic 
flight  rogimo.  An  udaptlvu-wul l  wind  tunnel1'1’ 
provides  for  uctlvo  control  of  tho  flow  in  the 
vicinity  of  tho  test-section  walls  in  order  to 
mlnlmizo  or  eliminate  tho  lnturforoncc  on  tho 
model,  the  distributions  of  tho  disturbance 
velocity  components  are  muanured  at  discrete  points 
along  control  surfaces,  or  Interfaces,  tn  the  flow 
fluid  near  the  walls,  but  away  from  the  model  and 
outside  the  boundary  layers  on  the  walls.  A  theo¬ 
retical  representation  for  tho  flow  exterior  to 
tlioso  control  surfaces,  Including  the  desired  uncon 
flned-flow  boundary  condition  that  nil  disturbances 
vanish  at  infinity,  Is  used  to  establish  tho  func¬ 
tional  relationships  which  must  bo  satisfied  at 
the  control  surfaces  by  tho  measured  disturbance 
velocities.  If  the  iuonsurod  velocities  do  not 
satisfy  tlioso  relationships,  an  Iterative  procedure 
provides  a  now  approximation  for  the  flow  field  at 
the  surfaces,  and  tho  wall  control  Is  readjusted 
until  tho  measured  quantities  satisfy  the  function¬ 
al  relationships  for  uncoil  fined  flow, 

Although  the  ultimate  application  Is  to  fully 
three-dimensional  flows,  most  of  the  adapt  1 vc -wa 1 1 
work  to  date  lias  boon  carried  out  In  two  dimensions 
Several  different  adapt lvo-wa 11  implementations 
have  boon  Investigated  experimentally  in  two  dimen¬ 
sions,  namely:  Impermeable,  flexible'  walls,11"11 
perforated  walls  with  constant  plenum  pressuro  but 
a  porosity  distribution  that  varies  In  tho  stream- 
wise  direction,1-  perforated  walls  with  constant 
plenum  pressuro  und  a  porosity  which  Is  uniform 
along  ea^h  wall,  but  can  vary  from  one  wall  to 
another1-  (Including  a  three-dimensional  applica¬ 
tion1'’),  slotted  walls  with  segmented  plenum 
chambers,1'1  and  perforated  walls  with  segmented 
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plenum  chambers.  In  the  latter  two  ventilated- 
w.ill  implementations,  the  pressure  in  each  of  the 
approximately  twenty  segmented  plenum  chambers  can 
bo  controlled  individually. 

Simulations  of  adaptive-wall  tunnels,  includ¬ 
ing  theoretical  representations  of  the  flow  within 
the  tunnel,  played  an  important  role  in  the  early 
development  of  the  concept.  Simulations  were 
carried  out  by  several  investigators  for  incom¬ 
pressible  flows, 4. 15  compressible,  subcritical 
flowslG-19  and  flows  which  are  supercritical  at  the 
tunnel  walls. d, 17-18  Not  only  did  these  studios 
demonstrate,  numerically,  that  the  overall  proce¬ 
dure  converged  to  unconfined  flow,  even  for  super¬ 
critical  walls,  hut  they  also  cast  a  great  deal  of 
light  on  the  nuturo  of  the  iterative  process, 
ospoclally  the  necessity  for  undorrelaxation  of  the 
itorutlons,  Furthermore,  the  accuracy  requlrod  in 
the  satisfaction  ot'  tho  functional  relationships 
for  unconfined  flow  was  examined.  Those  simula¬ 
tions  woro  highly  idoallzod,  however,  No  attempt 
was  made  to  model  tho  flow  in  tho  vicinity  of 
actual  wall  configurations.  Instead,  it  wns 
assumed  that  perfect  control  of  either  tho  stream- 
wise,  u  ,  or  normal,  -v  ,  disturbance  velocity 
componont  was  available  everywhere  along  tho 
control  stirfucos,  even  at  shock  wnvos,  Tho  only 
exception  was  u  study  by  Soars, l1-*  who  considered 
an  approximate  representation  of  linport'oct  wull 
control  and  examined  tho  orrors  in  tlte  convorgod 
solutions. 

A  two-dimensional  porforutud-wall ,  sogmontod- 
plonuin  configuration  luis  boon  investigated  oxtun- 
sivoly  at  Culspun  in  the  Ono-I'out  Wind  Tunnel  and 
tho  experimental  configuration  has  boon  described 
previously, 2-H  n  was  not  absolutely  necessary 
to  know  tho  dotulls  of  tho  flow  In  tho  vicinity  of 
tho  walls  in  order  to  achieve  unconfinod-f low  con¬ 
ditions.  That  Is,  If  tho  dasirad  disturbance 
velocities  at  tho  control  surfucos  could  bo  sot 
experimentally  at  each  stop  of  tho  itorutlvo  proce¬ 
dure,  tho  details  of  tho  flow  at  the  walls  wore  of 
llttlo  Interest.  However,  as  tho  oxportmontul 
investigation  proceeded, It  was  concluded  that 
oporutlonal  procedures  for  exorcising  wall  control, 
including  its  eventual  automation,  could  be  en¬ 
hanced  by  tho  development  of  more  realistic  simula¬ 
tions,  These  simulations  would  modol  tho  flow 
through  tho  wulls,  tho  Influence  of  that  flow  on 
tho  wull  boundary  layers  and,  consequently,  tho 
Intoructlor  with  the  flow  within  tho  control 
surfaces, 

Tho  simulation  methodology  developed  for  two- 
dlmonslonul  Incompressible  flows  is  described  In 
Section  II  with  some  details  given  in  Appendix  A. 
Results  for  tho  incomprossiblo-flow  approximation 
should  bo  representative  of  subsonic,  subcritical 
flows  In  general.  Results  for  the  baso-lino  simu¬ 
lation  of  a  61-blockngo  NAC.A  0012  airfoil  soctlon 
at  an  angle  of  attack,  a.  ,  of  4°  uro  prosontod  in 
Soctlon  III,  lnfluenco  functions  for  the  voloclty 
porturbutions  Introduced  In  tho  test  section  by 
changing  tho  pressure  in  ono  plenum  chamber  at  a 
time,  whllo  holding  tho  pressure  In  tho  other 
chamhors  constant,  are  presented  in  Section  IV  and 
are  compared  with  idoallzod  Influonco  functions 
which  are  derived  In  Appendix  1).  finally,  major 
conclusions  drawn  from  tho  simulations  are 
presented. 


11.  Methodology 
Tunnei  Configurator 

Tho  Cal  span  One-Foot  Wind  Tunnel  is  a  con¬ 
tinuous-flow,  closed-circuit  facility  that  operates 
at  Mach  numbers  from  about  0.5  to  0.95  at  a  unit 
Reynolds  number  of  2xl06  per  foot.  The  adaptive- 
wall  tost  section  is  two  dimensional,  with  perfo¬ 
rated  top  and  bottom  walls  of  22,5V  open-area  ratio. 
Thu  plenum  chambers  behind  the  perforated  walls 
have  bcon  divided  into  18  segments,  10  on  the  top 
and  8  on  the  bottom,  and  each  segment  is  connected 
to  a  pressure  and  a  suction  source  i' rough  indivi¬ 
dual  control  valves.  Tho  pressure  source  is  the 
tunnel  stilling  chamber,  and  the  suction  source 
is  an  auxiliary  compressor  discharging  into  tho 
tunnel  circuit  in  the  diffuser. 

Tho  model  und  test  soctlon  are  shown  schemati¬ 
cally  in  Fig.  1.  Tho  model  is  an  NACA  0012  airfoil 
section  with  a  6-inch  chord,  a  ,  (6V  solid  block¬ 
age)  and  is  situutud  midway  hotwoon  walls  that  uro 
located  at  ^  ■  *  h  »  *6.00  inches.  Tho  test- 
soction  instrumentation  at  the  control  surfaces 
consists  of  flow-angle  probos  und  static-pressure 
pipos.  Tho  control  surfaces  aro  locuted  outside 
tho  wall  boundary  layers,  namely  ut  ■  *J.9S  in. 
Tho  a  und  v  disturbance  velocity  components  aro 
determined  from  measurements  made  with  this  instru¬ 
mentation.^ 

Tho  simulation  proeudtii'.  that  is  developed 
boro  doos  not  provide  directly  for  tho  simulation 
of  tho  iterative  procoss  used  In  running  tho 
ilalspun  tunnel.  Thut  Is,  It  does  not  simulate  the 
plenum-pressure  adjustment  technique  that  Is 
used3"'’  to  sot  desired  (*-(*,  iy(>  distributions 
In  the  tunnel.  Such  simulations  could  ho  developed 
later,  If  desirod,  on  tho  basis  of  the  modol  des¬ 
cribed  hero. 

Calculation  Method 

Simulation  of  tho  flow  within  this  nduptlvo- 
wull  tost  section  roqulros  solutions  for  tho  flow 
over  the  airfoil  la  the  prosonco  of  tho  viscous 
effects  at  tho  perforated,  segmented,  plenum- 
pressuro-control led  walls.  Outside  of  the  wall 
boundary  luyors,  tho  flow  Is  assumed  to  ho  lnvisc.id. 
for  tho  lnvlscld  flow,  a  precede  u  ,s  required  for 
calculating  the  flow  ovor  a  model  with  proscribed 
distributions  of  tho  normal  voloclty  components  us 
boundary  conditions  at  tho  walls.  Two  tools  are 
noodod  to  treat  tho  wall  boundary  layers.  First, 
a  relationship  Is  required  hotwoon  tho  trunsplvu- 
tlon  velocity  at  tho  wall  and  the  pressure  drop 
across  the  wall.  Socond,  a  method  is  needed  for 
calculating  tho  characteristics  of  turbulont  bound¬ 
ary  layers  with  u  transpiration  voloclty  (cither 
Into  or  out  of  tho  walls)  which  Is  continuous  ovor 
ouch  plenum  segment,  hut  Is  likely  to  ho  disconti¬ 
nuous  at  junctions  between  segments, 

Invlscld-Plow  Modol 

A  computer  program  developed  originally  to 
predict  the  two-dtmonsionttl  flow  over  airfoils  in 
a  solid-wall  wind  tunnel  with  nonunlformly-shoarod 
froo  streums2(,‘2‘  hud  been  modifiod15  to  simulate 
the  invincld  flow  within  an  adaptivo-wall  wind 
tunnel  und  was  adapted  further  for  this  investiga¬ 
tion.  The  computer  program  can  accommodate  air¬ 
foils  with  arbitrary  thickness  and  cumber.  In  the 
program,  tho  airfoil  is  represented  by  a  vortex 
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distribution  along  its  actual  surface  while  the 
walls  are  represented  by  surface  source  distribu¬ 
tions.  These  vortex  and  source  distributions  are 
broken  up  into  s«all  segments,  each  of  which  has  a 
constant  strength.  Beyond  a  specified  distance 
from  the  airfoil,  both  upstream  and  downstream,  the 
source  distributions  at  tho  walls  are  assumed  to 
decay  inversely  with  t  approaches  ±  «  because 
the  source  strength  Is  closely  related  to  the  nor¬ 
mal  velocity  there.  The  volocity  components  in¬ 
duced  at  any  field  point  by  each  segment  are  given 
by  closed-form  algebraic  expressions.  Total  velo¬ 
city  components  due  to  the  walls  and  the  airfoil 
aro  then  found  by  summing  the  contributions  of  all 
segments.  The  strength  of  each  small  segment  Is 
found  by  satisfying  simultaneously  the  normal -flow 
boundary  conditions  v(X,ih)  at  the  wulls,  and 
tho  condition  of  no  normal  flow  through  tho  airfoil 
surfuco.  Thu  normal-flow  boundary  conditions  aro 
applied  at  tho  renter  of  ouch  segment.  For  N  vor- 
tox  and  sourco  segments,  N  algebraic  equations 
result  from  tho  boundary  conditions,  but  only  H- 1 
of  thorn  uro  linearly  indopendent,  us  shown  by 
von  Misus.22  \udit.lon  of  tho  Kutta-Joukowski  con¬ 
dition  completes  the  sot  of  equutlans,  which  is 
solved  using  u  smoothing  technique. 23-24 

Trunsplrod-Wull  Characteristics 

Tho  relationship  hotwoon  tho  transpiration 
vuloclty  at  tho  wall,  v«,  ,  (taken  as  positive  out 
of  tho  wull)  and  the  prossura  drop  uoross  the  wait 
was  dotorminod  empirically.  Data  wore  obtained  by 
Chow 2 5  for  tho  same  23.51  opon-urcu  rutlo  walls 
(1/10  la.  thick  plutu  with  i/16  in,  diumutor  holes 
on  1/S  in.  centers] ,  Those  dutu,  which  aro  only 
for  suet  lun,  wore  inousurod  at  Mach  numbers  from 
0.75  to  1.175.  Although  those  datu  wore  an  uvor- 
ugud  trw  over  a  long  run  of  pinto  with  constant 
plenum  pressure,  ,  and  uniform  prossura,  , 

In  tho  tost  suction,  It  was  assumed  thut  they  could 
ho  applied  In  tho  prosont  analysis  on  a  point-hy¬ 
po!  nt  basis,  Tho  data  of  Rof.  25  aro  corrolatod 
reasonably  wol 1  by 

(1°t  *  O.OOJi  +220.0  <V„/Um)'  (1) 

whore  <fa  and  M„  arc  tho  I’roe-stronm  dynamic  pros- 
suro  and  Mach  number.  This  result  had  to  ho  extend¬ 
ed  to  lower  values  of  und  to  Include  blowing. 
For  blowing,  Fig.  6,2  of  Hof.  26  Indicates  that  the 
data  could  reasonably  he  assumed  to  bo  untl-symmot- 
rlc  about,  vw«o,  Hof.  26  also  Lndleutos  that  at 
lower  Mach  numbers,  a  Frnndtl-lllauort  behavior  is 
expected.  Therefore,  a  I'rnndtl-Gluuort  scaling  was 
matched  to  the  empirical  variation  of  F.q.  (1)  and 
tho  anti -symmetric  character  wus  accounted  for  to 
obtain 

Cp„  *  <&■*>»)  /<}»  '  f(M)  [o,0OS -««.£>  |-^|(~^)] 
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Boundary-Layer  Modol 

The  most  important  characteristic  of  the  bound¬ 
ary  layer  for  the  present  application  is  its  inter¬ 
action  with  the  inviacid  flow.  Lighthill27  has 
described  several  ways  of  interpreting  the  displace¬ 
ment  thickness  and  its  effect  on  the  inviscid  flow. 
The  Interpretation  chosen  here  is  an  extension  to 
transpired  walls  of  what  iighthill  calls  the  method 
of  "equivalent  sources."  In  this  procedure,  the 
original  wall  surface  is  retained  and  the  affect  of 
the  boundary  layer  on  the  inviscid  flow  is  repre¬ 
sented  by  a  distribution  of  sources  on  this  surface, 
This  procedure  has  beoa  used  with  success  in  calcu¬ 
lations  of  tho  effect  of  tho  boundury  layer  on  the 
flow  over  airfoils  by  several  authors,  ax  in  Rofs. 

28  and  29,  for  example,  In  the  prosont  analysis 
thon,  the  effect  of  tho  boundary  luyer  on  tho 
inviscid  flow  is  written  as  an  inviscid  volocity 
distribution  14  (X)  (taken  as  posltivo  out  of  tho 
wall)  which,  when  applied  nt  the  actual  wull  sur- 
fac.0,  glvos  tho  proper  displacement  effect,  For 
trunxpirod  boundary  luyors,  v.  is  given  by 

ti  <*>  •  vw(*)  *d  /d x  (-*) 

whore  <V *  is  tho  boundary- layor  displacement  thick¬ 
ness,  and  U,  is  tho  total  Inviscid  stroumwiso 
volocity  component  outside  tho  boundary  layer.  In 
guttural,  <5*<x>  in  liq,  (4)  depends  on  u„ut),  lltIX> 
and  tho  Initial  conditions  for  tho  boundary- layor 
calculut ton, 

Since  tho  details  of  the  boundary- layer  profile 
aro  of  loss  concern  for  this  application,  Head's 
well-known  Integral  method**'1  was  adopted  for  the 
boundary- luyer  calculations .  An  outllno  of  the 
method  mid  representative  results  aro  given  In 
Appendix  A,  The  rositltlng  boundary- layer  equations 
aro  ittlugrntod  for  arbitrarily-proscribed  distribu¬ 
tions  of  Utt*i  and  Vw<jO  and  appropriate  ittltiul 
conditions. 

Tito  ostnbl  Isbmoat  of  tho  Invixcld-Flow  calcu¬ 
lation  method,  tho  wal 1 -character  1st lc  model,  and 
tho  boundnry-inyor  calculation  method  thus  sot  tho 
stage  for  development  of  a  simulation  procoduro 
for  approximating  uncon fined  flow  In  the  tunnel 
while  accounting  for  the  perforated  walls  with 
augmented  plenum  chambers.  First,  however,  I  clou  1 
pnlnt-by-polnt  wull  control  was  examined  from  the 
standpoint  of  transpired  turhulunt  boundary- layer 
behavior. 

goth  ideal  and  sogmontod-plomim  wull  control 
simulations  require  uso  of  the  invlscld-flow  pro¬ 
gram  to  calculate  the  unconflnad  flow  about  the 
model  ut  a  given  unglo  of  uttuck.  This  provides 
tho  unconflncd-f low  distributions  uis.ih)  and 
v  (X,  t  h  )  at  tho  walls  so  that  IV*>  *  l^  +  uf  X.th) 
and  bunco  ut  each  wall  are  known. 


whore 


f(M)  . 


5  O-mV* 

M 


M  t  (l)'1** 
M 


Ideal  Wall  Control 

Investigation  of  ideal  wall  control  wus  Ini¬ 
tiated  in  order  t.»  shod  light  011  how  the  wall 
boundary- layer  and  transpiration  volocity  churac- 
(;t)  toristlcs  must  Internet  to  provide  perfect  pnlnt- 
by-polnt  control  of  the  flow  within  the  tunnel. 
Calculations  with  this  procedure  then  serve  us 
guidelines  for  simulations  of  tho  imperfect  wull 
control  that,  exists  in  the  sogmented-plemim  config¬ 
uration. 


For  tho  incompressible-flow  calculations  of  this 
study,  f(M)  «0.5  wus  used. 
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The  unconfined-flow  distributions  u.(X.*hl 
and  v  (*,it i)  at  the  walls,  in  the  form  of  Ua(xi 
and  ,  ara  specified  on  a  point-by-point 

basis  in  order  to  calculate ‘ tha  ideal,  continuous 
wall  transpiration  velocity  distribution  v«  <*> 
which  would  reproduce  exactly  the  unconfined  flow. 
The  boundary- layer  problem  posed  by  this  ideal  wall 
control  Is  an  Inver so  one  because  Ut  and  v4  are 
proscribed  and  vw  must  be  found.  Because  of  the 
implicit  nature  of  thoir  dependence  on  vw  ,  tho 
governing  oquatlons  cannot  be  integrated  directly 
in  this  case,  and  an  iterative  solution  becomes 
necessary,  This  is  carried  out  at  each  step  of  the 
integration  us  follows.  A  value  of  vw  constant 
ovor  the  length  of  tho  stop  Is  assumed,  and  the 
aquations  uro  integrated  for  tho  stop.  Eq,  (4)  is 
ovuluatod  noxt  to  determine  tho  calculated  v«  which 
is  eompnrod  with  tho  desired  value.  If  they  do  not 
agree  to  a  specified  ncourucy,  tho  value  of  vw  Is 
roudjustod  and  tho  process  is  repeated  itorativoly 
until  tho  calculated  vt  agrees  with  the  proscribed 
vnluo.  The  v„iri  distributions  which  result  from 
this  calculation  provldo  a  useful  roforonco  for 
beginning  moro  realistic  sogmontod-plonum  calcula¬ 
tions,  as  shown  noxt. 

Sogmontod-Plonum  Wall  Control 

Thu  modeling  of  sogmontod-plonum  wall  control 
also  rollos  on  tiiu  unconflnod- f low  distributions 

and  V,  (X)  .  The  hourt  of  tho  procedure 
Is  thun  tho  dotormlnution  of  +>„  ,  tha  plonum 
pressure  In  thu  n  ™  plenum  chamber,  This  Is  accom¬ 
plished  Iteratively  as  follows.  At  tho  contur  of 
each  plonum  chamber,  x„  ,  tho  invlscid-flow  pres¬ 
sure,  Is  known,  and  since  un  Initial 

estiumto  fur  vw(x„)  cun  bo  obtuLnod  from  tho  ldonl- 
wall  calculation,  un  initial  guoss  for  cun  bo 
found  from  Uqs.  (2)  und  (3),  Kith  and  (X ) 
thun  known,  tho  corresponding  Initial  approximation 
to  11m,  (X)  is  found  ovor  tho  length  of  ouch  plenum 
clumber  by  !!q.s.  (21  und  (3).  In  general,  this 
procedure  loads  to  discontinuities  In  vw  <*>  ut 
tho  boundaries  butwoen  plenum  chambers,  Thu  hou:  <1- 
nry-luyor  calculation  is  thon  made  ovor  each  plenum 
clumber  In  turn  on  tho  basis  of  tho  proscribed  In¬ 
put  O,  tr„(X),  und  tho  Initial  conditions.  Tho 
offoctlvo  Invlscid  velocity  distribution  at  the 
walls,  v,  (X) ,  U  found  from  tho  computed  results 
and  llq.  (4).  Tho  Integral  of  v,  <.%)  ovor  each 
plonum  chamber  is  thon  compared  with  the  lntogrul 
of  thu  unconflnod" flow  distribution  V  (X,  th) 
ovor  tho  same  plenum.  (It  should  bo  noted  that  In 
tho  dlfforont  coordinate  systems  used  for  thu  In- 
viscid-flow  tmd  boundary- layer  calculations, 
v  <  x,h  1  «  -v,  (t )  und  v  <*,-h)  ■  vt(x)  .) 
Tho  ■?„  nro  thon  revised  and  tno  boundary- layer 
calculation  Is  ropoatod  until  tho  integrals  of 
vt  (t )  und  v(x,ih)ovor  ouch  plenum  ugroo.  Satis¬ 
faction  of  this  condition,  locally  at  each  plenum, 
insures  that  tho  lntogrul  of  the  source  olomants 
ovor  ouch  wall  Is  zero  us  is  tho  case  for  fully 
unconflnod  flow.  This  Is  equivalent  to  tho  fact 
that  tho  drag  on  thu  airfoil  Is  zoto  in  unconflnod, 
invlscid  flow. 

Onco  tho  VfcOO  distributions  nt  both  walls 
have  boon  found  In  this  way,  tlioy  arc  used  as  bound¬ 
ary  conditions  In  tho  Invlscid- flow  program  to  cal¬ 
culate  tho  flow  ovor  tho  airfoil  In  tho  presence  of 
tho  sogmontod-plonum,  perforntod-wnll  constraints. 
Ilcsults  found  by  such  a  procedure  can  he  usod  to 
Invest Ignto  tho  magnitude  of  tho  orrors  Introduced 
by  having  segmented- plonum,  lmporfoct  wall  control 


instead  of  ideal,  perfect  wall  control;  this  is 
discussed  in  the  next  section. 

III.  Baso-Llne  Simulation  Results 

The  example  chosen  for  the  base-line  simulation 
is  the  NACA  0012  airfoil  of  the  experiments  ut  an 
angie-of-attack,  «.  ,  of  4*.  This  choice  gives  a 
good  balance  between  tho  lift  and  thickness  effects 
on  the  disturbance  velocity  field  at  the  control 
surfaces  and  walls.  Wall  and  control -surface  loca¬ 
tions  at  1.0  and  0,655  chords  above  and  below  the 
model  correspond  to  the  experimental  configuration. 
Within  this  geometrical  Input,  both  ideal  and  seg- 
monted-plonum  wall  control  were  examined. 

for  ideal  wall  control,  thu  calculated  distri¬ 
bution  of  and  the  prescribed,  unconfi nod- 

flow  distribution  «/*<*>  are  presented  in  fig.  2 
for  tho  uppor  wall.  Note  that  upstroum  of  the 
origin  thoro  is  suction  und  a  favorable  pressure 
gradient  so  v4(x>  follows  vw(X)  closely.  Down¬ 
stream  of  tho  origin  thoro  is  blowing  und  an  ad¬ 
verse  pressure  gradient ,  which  givo  rise  to  n  ratio 
of  Vg  to  vw  of  two  or  grouter,  i.o.,  an  amplifi¬ 
cation  of  tho  wall  velocity  when  blowing  Is  present. 
Such  behavior  is  typical,  nnd  tho  rousons  fur  it 
nro  described  nt  groutor  length  In  Appendix  A. 

Tho  results  of  fig.  1  weru  usod  to  inltluto 
the  Iterative  calculation  of  thu  struumwisa  distri¬ 
butions  of  V„(*l  for  tho  sogmontod-plonum  model, 
as  described  in  Section  11.  Tho  rosult  for  tho 
uppor  wall  Is  shown  us  the  solid  curve  in  fig.  3; 
results  for  tho  lower  wall  dlspluy  u  similar 
behavior.  The  tick  marks  on  tho  ubsclssa  of  fig.  3 
correspond  to  tho  Junctions  between  plonum  chambers 
as  shown  In  fig.  I.  for  comparison,  tho  ideal 
polnt-by-polnt  distribution  of  ,vw(X>  from  fig,  2 
Is  shown  us  tho  dashed  curve  on  fig.  3.  Tho  nddl- 
tlonnl  curves  und  data  points  on  this  and  the  noxt 
few  figures  will  he  discussed  In  Soctlon  IV.  Tho 
stroumwiso  distributions  of  v„<%)  thut  rosult 
from  tho  matching  procedure  have  very  largo  varia¬ 
tions  ovor  each  plonum  chamber,  with  discontinuities 
at  tho  Jin.  t.ions  botwoon  plonum  chambers, 

Thu  discontinuous  distributions  yield 

correspondingly  largo  variations  In  va  from 
the  boundary- layer  calculations,  particularly  as 
amplified  when  blowing  is  present.  Tho  distribu¬ 
tion  of  v4<<>  for  the  upper  wull  Is  given  us  the 
solid  curve  In  fig.  4,  along  with  tho  dashod  curve, 
which  is  the  unconflnod- flow  distribution  vCx,b) 
from  fig.  2.  Tho  largo  variations  in  Vt<x>nnd 
tho  associated  di  scontimiltios  at  tho  plonum  Junc¬ 
tions  arc  certainly  not  roullstlc  physically.  The 
response  of  tho  boundary  layers  at  tho  junctions 
would  not  havo  such  un  oxtromo  character.  Never¬ 
theless,  those  details  do  not  have  a  significant 
Influonco  on  thu  resultant  flow  conditions  at  tho 
control  surfaces  or  model,  us  will  bo  soon  below, 
so  this  representation  was  retained. 

Tho  scgmontod-plcnum  v4<x>  distribution  of 
fig.  4,  and  Its  counterpart  ut  tho  lower  wall, 
wore  usod  as  boundary  conditions  In  the  calculation 
of  tho  Invlscid  flow  about  tho  airfoil,  Rosults 
for  the  stroumwiso  and  normal  velocity  components 
at  tho  uppor  and  lower  control  surface*  uro  present¬ 
ed  in  figs,  5  und  6  by  tho  crossos.  The  unconfined- 
flow  distributions  nro  given  by  the  dashed  curvo 
while  the  rosults  given  by  tho  circles  will  he  dis¬ 
cussed  in  Section  IV. 
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Examination  of  Figs.  5  and  6  indicates  that 
tho  agreomont  between  the  segmented-plenum  and 
unconfined- flow  distributions  is  generally  excel¬ 
lent  where  suction  is  present  at  the  walls,  namely 
upstream  of  the  origin  at  the  upper  control  surface 
and  downstream  of  the  origin  at  the  lower  control 
surface.  The  discrepancies  observed  in  u  at  the 
upstream  and  downstream  limits  of  the  test  section 
result  from  the  truncation  of  v*  it)  from  the  dis¬ 
tributions  shown  in  Fig.  4  to  zero  upstream  and 
downstream  of  tho  controlled  part  of  the  test- 
section  walls.  Repeating  the  inviscid  flow  calcu¬ 
lation  with  Vet*)  oqual  to  tho  unconfinod-flow 
distributions  boyond  the  controlled  sections  result¬ 
ed  in  excellent  agreement  of  a  at  the  upstream  and 
downstream  limits  of  tho  test  soction, 

Overall,  tho  results  indicate  that  tho  finite 
slzo  of  tho  plonum  chambers  und  th«  attendant  largo 
variations  In  v«  cx>  from  modeling  tho  flow  ut  tho 
walls  do  not  prevent  tho  realization  of  a  good 
upproximutlon  to  tho  unconfinod-flow  distributions 
for  tho so  low- speed  flow  conditions.  This  is  con¬ 
sistent  with  the  onrllor  experimental  demonstrations 
of  tho  adaptive-wall  concept  for  flows  which  aro 
supercritical  at  the  model,  but  subcvitical  at  tho 
control  surfaces  and  walls, ^i4  in  thoso  casos,  the 
desired  flows  could  ho  sot  up  experimentally  with¬ 
out  largo  fluctuations  In  tho  measured  stroumwise 
volocity  components  from  point  to  point. 

Thu  dl strlbut Ions  of  the  pressure  coefficient, 

,  on  tho  uirfoll  surface  lire  given  In  Fig.  7, 
with  tho  dashed  linos  indicating  unconflnod  flow 
and  tho  crosses  tho  results  of  tho  sogmonted-plonum 
model,  Tho  excellent  agreement  shown  Is  confirmed 
by  integrating  the  distributions  to  obtain 
normal  force  coefficients  of  0.490  for  unconflnod 
flow  and  (),4Hd  for  tho  stpnontod-plonum  model,  a 
dlffuronco  of  only  1,6'!.,  Tho  remarkable  agreement 
In  both  the  airfoil  Cp  distributions  and  the  dis¬ 
turbance  volocity  distributions  at  tho  control  sur¬ 
faces,  despite  the  unrealistic  nature  of  tho  v€  <x> 
distributions,  results  from  tho  elliptic  nature 
of  the1  governing  equations  of  motion  for  tho  flow 
within  the  tunnel,  Furthermore,  this  behavior  is 
roproxuntnt 1 vo  of  tho  uctuul  tunnel  flow. 

Wo  huliuvo  that,  thoso  simulations,  dosplto  tho 
exaggerations  in  tho  variations  of  vw  <*)  and 
vt<x)  over  tho  individual  plonum  chambers,  espe¬ 
cially  ut  the  plonum  junctions,  model  tho  essential 
features  of  tho  flow  nour  the  walls.  Wo  holiovo 
further  that  the  qualitative  features  of  those 
huso- lino  simulations,  although  they  aro  strictly 
for  incompressible  flow  only,  can  bo  generalized 
to  compressible,  suhcrltlcul  flow  because  there  are 
no  fundamentu 1  differences  In  tho  flow  phonomonu 
Involved.  However,  onco  tho  flow  becomes  super¬ 
critical  ut  tho  wulls,  dlfforont  phonomonu  are  In¬ 
volved  and  modeling  of  tho  shook-wuvo/boundary-lnyor 
Interaction  becomes  necossury. 

Tho  existing  simulation  methodology  can  serve 
us  the  basis  for  further  investigations  of  tho 
behavior  of  udaptlvo-wull  wLnd  tunnels  with  this 
wall  configuration.  For  example,  it  can  bo  used  to 
examine  tho  effects  jf  plonum  chamber  number  und 
size  and  of  the  porosity  distribution  of  the  perfo¬ 
rated  wulls.  Also,  it  forms  a  frumowork  for  t.ho 
development  of  more  systematic  procedures  for 
adjusting  tho  tunnel  flows  that  uro  required  ut 
ouch  stop  of  the  overall  uduptive-wul 1  itorutivo 
process.-’-'1  Thoso  adjustment  procedures  would 


provide  the  logic  for  the  ultimate  automation  of 
adaptive-wall  wind  tunnels.  A  first  step  toward 
logical  procedures  for  tunnel  adjustment  is  des¬ 
cribed  in  the  next  section. 


IV.  Influence  Functions  Due  to 
Plenum  Pressure  Change 

As  a  first  step  toward  the  development  of 
automated  wall-adjustment  procedures,  u  study  was 
made  of  the  effect  on  tho  flow  produced  by  chunging 
the  prossure,  +n  ,  in  one  plonum  chamber  ut  a 
time,  while  holding  the  pressure  in  the  othor 
chumbors  constant.  Tho  oxumplo  doscrlbod  in  Soction 
III  was  chosen  for  those  calculations.  Numerous 
casos  were  investigated,  including  plonum  chambers 
at  various  stroumwise  locutions  in  tho  tunnol  with 
suction  and  with  blowing. 

In  one  pair  of  examples,  the  prossure  in 
plonum  6  (sec  Fig.  1)  was  changed  by  amounts  such 
that  v^ix^)  /  was  changed  by  *0.0100  (und  so 
^  -0.051  and  0.053,  rospoctivoly) ,  where 
x„  is  tho  *  -coordinate  of  tho  center  of  plonum 
m  .  The  ovorull  changes  in  v„cx)  over  plenum  6 
urn  shown  In  Fig.  3.  The  v«  u>  distributions 
which  result  from  the  boundary -layer  calculation 
nro  shown  on  Fig,  4.  It  Is  romurkublu  that  those 
resulting  v,(t)  distributions  uro  changed  signi¬ 
ficantly  along  the  longtli  of  that  plenum  chamber, 
but  negligibly  at  all  other  chambers,  even  those 
Immediately  downstream.  This  behavior  occurs  for 
nil  oxumples  consldc-rod,  whether  the  basic  condi¬ 
tion  Is  suction  er  blowing  at  the  plenum  chamber 
where  tho  prossure  change  Is  made.  This  does  not 
Imply  that  tho  boundary  layer  Is  unaffected  down¬ 
stream.  Many  of  Its  properties,  o.g.,  6*  und 

<9,  change  appreciably  down  at  ream  |  however, 
diUt< S‘)  /  dt  and  hence  v«(x>  from  liq.  (4)  are 
not  affected  significantly  there. 

Results  from  the  Invlscld-flow  calculation 
for  the  Increased  suction  example  of  Fig.  4 
(  A  trvv<  xu)  /  Ua  »  -i).  01 00  or  «  0,05.3)  uro  shown 

ns  tho  circles  In  Mgs.  5  and  o,  Whim  compared 
with  tho  base-lino  simulation  represented  by  the 
crosses,  It  Is  observed  In  Fig,  5  that  u-  hits  In¬ 
creased  upstream  and  decreased  downstream  of  plenum 
(i  ut  both  control  surfaces  with  u  crossover  nt 
plenum  6  Itself,  This  huhnvlor  Is  consistent  with 
the  presence  of  a  sink  nt  the  wall  of  a  confined 
channel,  such  as  exists  here,  An  Increase  In  suc¬ 
tion  (or  ii  docroaso  In  blowing)  at  one  plonum 
chamber  gives  rise  to  an  Increased  sink  strength 
nt  that  wall  location.  Conversely,  other  results 
not  prosontod  hero  show  that  a  decrease  In  suction 
(or  an  increase  in  blowing!  gives  rise  to  an  In¬ 
creased  source  strength  and  tho  differences  In  tho 
u-  distributions  will  bo  of  opposite  sign.  This 
source  (or  sink)  offoct  In  un  otherwise  confined 
channel  Is  felt  oven  to  large  distances  upstream 
and  downstream,  as  observed  in  Fig.  5. 

Tho  sink  offoct  described  here  was  observed 
in  an  experiment  performed  In  I  lx-  Cal  span  Ono- 
Foot  Adnptlve-Wttl 1  Tunnel  before  this  analysis  was 
undertaken.  Tho  suction  wus  Inc  run. soil  simultaneous¬ 
ly  ut  both  plonum  chumbors  5  and  15,  Tho  increased 
suction  increased  tho  flow  velocity  at  that  station 
far  upstream  whore  M„  is  measured  and  set,  In 
order  to  maintain  the  doslrod  M„  ,  compensating 
adjustments  wore  made  In  tho  tunnel  drlvo,  effec¬ 
tively  removing  tho  upstroum  perturbation. 


Therefore,  during  experimental  iterations  at  a 
constant  M„  ,  one  can  expect  the  primary  changes 
in  u  to  occur  at,  and  downstream  of,  the  plenum 
being  adjusted. 

On  the  other  hund,  the  chunges  in  the  v  dis¬ 
tributions  at  the  control  surfaces,  as  shown  in 
Pig,  6,  occur  only  locally  in  the  immediate  vicin¬ 
ity  of  the  plenum  chamber  at  which  the  control  is 
changed.  There  is  only  a  slight  effect  at  the 
same  stroamwlse  position  on  the  opposite  control 
surface  and  negligible  effect  upstream  and  down- 
stroam  at  both  control  surfaces.  This  localized 
behavior  also  was  observed  in  the  Ono-Foot  Tunnel 
experiments. 

The  rosultn  of  all  the  calculations  in  which 
thu  pressure  in  u  single  plenum  chamber  was  changed, 
while  holding  the  remainder  at  their  original 
pressures,  wore  analyzed  und  cast  in  terms  of 
lnfluonco  functions.  Thu  notution  Au.<-xi  will 
rotor  to  the  difforonco  given  by  tho  streamwiso 
velocity  porturbution  uftur  a  plenum  adjustment 
has  heon  mude  minus  its  vnluo  in  tho  base-line 
simulation;  un  analogous  definition  is  made  for 
A  v  tx  >  .  If  fa.  i  '  aa('»)  is  tho  differ¬ 

ence  in  this  calculated  perturbation  aver  tho 
length  of  tho  test  section,  then  a  naturui  form  for 
expressing  tho  Influence  functions  Is 

Iu<£)  -  [au<S:)  -  nut-*)]  /6a  (S) 

lvt  it)  «  -\Aviz)/6a\  ('>) 

whore  X  «  (  x-x„)/Jh  1  a  a  coordinate  normalized  by 
in  ,  the  length  of  plenum  m.  .  Tho  normal  l  tut  ton 
in  lips.  I5j  and  (0)  reflects  tho  oxporlmoatal  situa¬ 
tion  In  which  the  tunnel  operator  removes  any  A u. 
that  Is  introduced  fur  upstream  so  os  to  maintain 
a  constant  Mw  i  hence  "  U  and  Iu«»)  •  1. 

The  form  of  Iv  ,  Including  the  sign,  Is  chosen  so 
that  pressure  changes  at  both  upper  and  lower 
plonum  chambers  uve  cons  1  stout  with  tho  results  of 
Appendix  I) .  (Sec  tho  discussion  loading  up  to 
h-.|s.  (11-19)  and  lll-:o).'l 

Tho  lnfluonco  functions  I  u,fc  and  ,  where 

flu-  subscript  a.  denotes  the  control  surface 
adjacent  to  pieman  n  ,  are  presented  In  figs.  H  to 
II  as  functions  of  x  ,  where  the  lines  nt 
x  =>  *0.5  Uemnreato  the  ends  of  nleiitun  n  .  In 
Figs.  8  and  !>,  tho  data  from  all  tho  plenum  cluiiiihers 
considered  In  tho  calculations  are  plottod,  namely 
plenum  chambers  d,  S  and  fi  with  incrousod  suction, 
plonum  ()  with  decreased  suction,  and  plenum  ch.Uiiburs 
II,  1.1  and  In  with  Increased  blowing.  For  compari¬ 
son,  Idonllzod  Influence  functions,  j und 

lr,  i  uro  plottod  as  solid  curves.  These  Ideal¬ 
izations  itro  derived  la  Appendix  It  using  linearized 
theory  without  a  model  present  and  assume  a  uniform 
Ai>t(iO  between  i  •  *0. 5.  Figs.  10  and  II  are 
plotted  to  an  on larged  x  scale  to  Illustrate  how 
t.lio  details  of  the  changes  in  the  vt  distribution 
in  lUg.  4  nffoct  tho  lnfluonco  functions  in  tho 
Immediate  vicinity  of  plonum  (>.  .Since  the  lnc.ro- 
monts  lit  v.  In  that  vicinity  arc  larger  at  tho 
downstream  end  of  tho  plenum  than  at  tho  upstream 
end,  and  In  Figs.  10  und  11  aru  skewed 

us  shown.  Further  uwuy  from  thu  plonum  tho  Idoul- 
1 zod  Influence  functions,  which  ussumo  u  constant 
along  tho  plenum,  provldo  a  satisfactory 
approach  to  tho  asymptotic  conditions. 


The  Influence  functions  1^  and  ,  where 

the  subscript  o  denotes  the  control  surface  at 
the  opposite  side  of  the  tunnel  from  plenum  n  , 
are  presented  in  Figs.  12  and  11.  Again,  the  de¬ 
tails  of  the  boundary- layer  behavior  at  plenum  n 
primarily  affect  the  distribution  only  in  the  imme¬ 
diate  vicinity  of  that,  plenum  while  the  behavior 
away  from  this  region  is  approximated  joasonahly 
well  by  the  idealized  representation,  ltt(  and 

.The  fact  that  the  idealized  distributions  i„, 
and  1„  derived  for  linearizod,  empty-tunnel  flows 
are  such  a  good  representation  of  the  distributions 
computed  with  the  model  present,  1^  und  Iv  ,  is 
quite  romurkable.  As  can  be  seen  from  the  deriva¬ 
tion  of  the  linourizod  model,  the  general  behavior 
is  a  function  of  plenum  length  und  the  distnneo  of 
tho  control  surfaces  from  tho  wulls.  It.  is  expect¬ 
ed  that  other  sogmontod-plonum  configurations,  o.g., 
thoso  with  slotted  walls  or  those  with  porforuted 
wulls  of  different  opon-areu  ratios  or  aluntod 
holes,  should  exhibit  tho  sumo  general  characteris¬ 
tics, 

The  relationship  between  the  change  in  pressure, 
say,  in  plonum  n  and  the  total  change  in 
tho  u.  component,  <Ju,  ,  will  depend  on  thu  wall 
geometry,  of  eourso.  This  relationship  depends  on 
the  details  of  thu  bound.,  •*>  layer  and  tho  wall 
characteristics  und  for  the  dulspnn  One-Font  Tunnel 
wulls  Is  Implicitly  repress  n*-ed  by  tho  form  of  llqs. 
(15  to  (4).  For  example,  In  , ho  cases  presented  In 
Figs.  1  and  4,  the  Increment  zlC^  »  -(1.0310  results 
In  du -  0,0009  and  tho  Increment  iC,,k  «  0.0810 
results  In  t5u/U„"  -0,0017,  whore  these  cases 
represent  Increments  In  A vw(xu>/U<,of  0.0100  and 
-0,0100,  respectively. 

There  are  two  Important  next  stops  in  the 
development  of  this  simulation  and  Influonco- 
functlon  technique.  The  first  would  ho  the  exten¬ 
sion  of  the  analysis  to  compressible  flows,  As 
mentioned  In  Section  111,  no  significant  differences 
lire  expected  tint  1 1  supercritical -flow  condition!! 
are  reached.  Thou  ll  Is  expected  that  T,t  and  Sv 
would  depart  significant lv  from  the  Idealized  sub¬ 
sonic  behavior.  The  other  stop  would  ho  tho 
experimental  verification  of  tho  In  flounce  -  fund  ion 
bulmvlor  by  measuring  tho  npproprli.  o  quantities. 
This  would  he  straightforward  hut  was  beyond  the 
scope  of  tho  experimental  programs  porformod  thus 
far.  Ilodiiputi,  et  til.l'l  have  measured  Isanti  Zv, 
for  empty-tunnel  conditions  In  their  sogmontod- 
pleiutm,  slottod-wnll  test  section,  Their  results 
are  generally  of  tho  form  predicted  hero. 

Once  tho  validity  of  the  predicted  Influence- 
function  behavior  has  boon  established  explicitly, 
the  procedures  could  ho  used  to  investigate  auto¬ 
mated  control  procedures.  For  example,  in  the 
dal  span  experiments,  the  u.  component  is  set  along 
the  control  surfaces  so  that  in*  and  la.  are  of 
principal  importance.  In  order  to  effort  a  change 
in  u  at  oacii  control  point,  tho  sum  of  the  chrngus 
Induced  nt.  that  point  duo  to  pressure  changes  AC^ 
lu  each  of  tho  i\  plenum  chambers  must  bo  found. 
Since  similar  changes  are  desired  at  ull  the 
centre)  points,  a  system  of  equations  can  bo 
derived  relating  the  prescribed  changes  in  u.  to 
tho  plonum  prossuro  changes.  This  system  of  equa¬ 
tions  must  he  solved  for  thu  A Cp  that  are  nocos- 
sary  to  achieve  the  desired  flow."  Since  tho 
influence  functions  of  F.qs,  (5)  and  (6)  are  defined 
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in  terms  of  £1*  ,  this  quantity  must  be  related  to 
the  pressure  change  for  each  plenum.  As 

discussed,  this  relationship  is  a  highly  nonlinear 
one.  Use  of  both  the  idealized  and  tho  more  exact 
relationships  could  be  investigated  by  simulations 
to  establish  the  accuracy  to  which  the  influence 
functions  must  be  represented  in  the  vicinity  of 
each  control  plenum  to  iterate  automatical ly  and 
quickly  to  unconfined  flow. 


V.  Concluding  Remarks 


Tho  following  are  viewed  as  the  principal 
results  and  conclusions  drawn  from  the  present 
Investigation. 


To  the  authors'  knowledge,  the  base-line  simu¬ 
lation  represents  tho  first  realistic  numerical 
demonstration  that  an  adaptive-wall  wind  tunnel 
with  imporfoct  wall  control  can  nevertheless 
achiovo  what  is,  for  most  practical  purposos, 
uncoil  1‘lnod  flow  about  the  modei  (soo  Fig,  7).  Tills 
is  uccompl l shod  in  spite  of  tho  fact  that  tho 
modeling  of  the  wall  transpiration  characteristics 
and  houtidury  layors  was  relatively  crude,  and  if 
anything  oxacorbntod  tho  effects  of  imperfect 
control . 


Thu  local  boundary- layer  displacement  effect., 
as  embodied  by  Ve<*Hin  liq.  (4),  Is  influenced  pri¬ 
marily  by  pressure  chungos  In  the  plenum  chambors 
In  tho  linmodlato  vicinity,  and  Is  somewhat  depend¬ 
ent  on  tho  detailed  distributions  of  transpirution 
voloocy  through  the  walls  (soo  Pig.  4). 


When  appropriately  defined,  tho  Influence 
functions,  1*.  und  Iv,  appear  to  fall  on  a  more  or 
loss  universal  curve  lndopondont  of  plenum  location 
und  whether  suction  or  blowing  is  being  applied 
(soo  Pigs.  8-13).  However,  tho  relationship  botwoon 
AC*„  and  tho  normalization  factor  (liqs,  (5)  und  (6)) 
will  depend  on  the  particular  tunnel  configuration. 


Tho  influence  functions  well  away  from  tho 
immediate  vicinity  of  u  given  plenum  appear  to 
depend  only  on  tho  integrated  normal  velocity 
and  muy  bo  prodicted  adequately  by  a  simplified 
theory. 


Though  tho  presont  simulations  arc  confined  t.o 
low-speed  flow  for  simplicity,  tho  extension  to 
comprossiblo  subsonic  spoods  is  not  oxpoctod  to 
ultor  thoso  conclusions  qualitatively,  provided 
tho  flow  remains  subcrltlcal.  Por  supercritical 
wull  conditions,  which  muy  prove  of  groat  practical 
significance,  a  modol  of  sboek-wavo/boundary- layer 
interaction  would  have  to  bo  addodj  tho  attendant 
possibility  of  separated  flow  at  tho  wull  und  the 
sensitivity  of  shock-wave  position  to  details  of 
tho  flow  at  tho  walls  may  load  to  qualitatively  now 
phenomena.  Once  this  quostlon  is  answered,  tho 
application  of  the  model  to  simulating  tho  untual 
Itorativo  wall-control  procoss^'d  at  high  subsonic 
spoods  should  bo  possible. 
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Append  1.x  A 

Boundary- Layer  Model  and  Amplification  Effect 

l'or  thu  purpose  of  studying  tho  effects  of  tho 
boundary  layer  on  the  tnviscid  interior  flow,  it  is 
tho  integral  properties  of  the  layer,  us  opposed  to 
tho  dotailod  profiles  across  It,  which  uro  of  pri¬ 
mary  Interest.  Accordingly,  an  intogral  (as  opposed 
to  flnitc-difforcnco)  method  was  chosun  us  tho  most 
efficient  moans  to  model  tho  boundary  layor.  The 
particular  method  used  Is  that  first  proposed  by 
I  loud ,  which  had  been  applied  successfully  by 
others  In  predict Ingtho  behavior  of  transpired 
turbulont  layers. 51-32 

Basically,  Hoad's  method  consists  of  simulta¬ 
neously  solving  three  equations.  The  first  is  his 
so-called  entrainment  equation, 

±  jL  (U,6H,)  •  (A-l) 

<4 

where  is  a  modified  shape  factor  defined  by 
H,  ■  C and  U*  Is  tho  local  stroumwlso 
velocity  at  tho  layer's  edge,  Hero  6,6*  ,  and 

9  rotor  respectively  to  tho  boundary- layor,  dis¬ 
placement,  und  momentum  thicknesses.  H,  is  rela¬ 
ted  to  tho  more  conventional  shape  factor,  H  «  Sk/6, 
by  H,  •  G(W>  .  The  functions  8  (H,l  und 
Gr IH)  wore  empirically  derived  from  existing 
data.  Tho  rationale  hohind  Eq.  (A-l)  can  ho  found 
in  Refs,  30-32. 


The  second  equation  Is: 


i£  -  +  - 

d  X  Ut  dx 
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which  is  simply  tho  well  known  von  Kurmun  integral 
momontum  equation  derived  in  most  boundary  layer 
texts.  The  only  now  quantity  introduced  is  Cf  , 
tho  local  skin-friction  coefficient. 
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Head's  original  applications  ware  to  solid- 
wall  boundary  layers,  for  which  he  used  the  Ludwieg- 
TillMnn  correlation*3!  to  predict  Cf  .  For  the 
present  application  this  procedure  was  Modified  as 
follows:  first,  the  Ludwieg-Tillmann  correlation 
was  used  to  compute  Cf  ,  the  akin-friction  coeffi¬ 
cient  that  would  obtain* if  no  transpiration  were 
present , 


Cft  -  (0.E4U)  10 


-0.4TJW 


where  fo*  is  the  local  momentum-thickness  Reynolds 
number.  This  value  was  then  corrected  for  the 
effects  of  blowing  or  suction  via  the  empirical 
correction  suggested  by  Simpson,  Moffat,  and  Kays:34 


c,  ri«.  | ;+  5|  i 0,7 

**,  L  D  J 


(A-d) 


where  £1*  2vm/CfUt.  Due  to  the  implicit  depend 
ence  of  0  on  Cf  and  the  transcendental  nature 
of  Eq.  (A-4),  it  had  to  be  solved  using  Newton- 
Raphson  iteration,  but  the  convergence  is  quite 
rapid. 


The  above  system  of  throe  equations  (two  dif¬ 
ferential  ,  one  algebraic)  wus  solved  for  the  three 
unknowns  6  ,  H  and  Cf  using  a  Runge-Kutta 
procedure.  Required  inputs  aro  the  distributions 
ll4  and  -u’w<jo  ,  and  initial  conditions  on 
8  and  H  at  the  starting  point. 


To  vulidate  the  boundary- layor  program,  its 
predictions  were  comparod  with  oxpot'irnontal  data  of 
other  investigators  for  several  representative 
flows,  including  various  stroumwiso  pressure  gra¬ 
dients  with  and  without  transpiration.  An  oxamplo 
which  is  of  particular  interest  to  the  present 
application  is  un  experimental  configuration  inves¬ 
tigated  by  McQuuid.3*  In  this  cusc,  there  was 
uniform  blowing  of  constant  strength  up  to 
X  ■  17.5  in.;  downstream  of  this  point  nr*  wus 
nominally  ;oro.  The  equations  were  Integrated  up 
to  X  «  17.5  in.  using  tho  oxporimontal  vuluos  of 
8  and  H  at  X  -  11.5  in.  as  initial  conditions 
and  ■  0.0(135.  Tho  calculated  values  of 

8  and  H  at  x  »  17.5  in.  woro  then  used  as 
Initial  conditions  for  further  integration  down¬ 
stream  with  V+,  /U-“  0.  Results  of  tho  calcula¬ 
tion  are  s  lown  in  Fig.  14.  Tho  predicted  develop¬ 
ment  of  6  *  ,  8  and  H  from  their  initial  vuluos 
agrees  well  with  tho  data.  In  particular,  tho 
culculutod  distributions  display  discontinuities  in 
their  slopes  at  tho  discontinuity  In  blowing  velo¬ 
city,  and  his  behavior  reflects  tho  observed  varia¬ 
tion  there  reasonably  woll.  This  and  other  compar¬ 
isons  gave  confidence  in  tho  ability  of  thu  model 
to  predict  the  intogrul  properties  of  the  boundary 
layor  aduq  lutoly. 


The  f  rst  application  of  tho  boundary-layer 
model  specifically  to  tho  Culspun  One-Foot  Tunnel 
configuration  shown  in  Fig.  1  was  dono  to  see 
whether  It  confirmed  a  so-called  "boundary- layor 
amplification"  effect  that  hud  boon  noted  in  adjust¬ 
ing  the  plena  to  achieve  a  dosired  flow  condition. 

In  early  experiments  with  tho  model  present  in  the 
tunnel, 4  it  was  observed  in  many  cusos  that  the 
normal  velocity  measured  by  a  probe  l.S  in.  from 
the  wall  was  considerably  gruater  in  magnitude 
than  Vw  ,  tho  velocity  at  the  wull  itself. 

Although  part  of  this  amplification  was  due  to  the 
increase  ovor  that  diatance  Ln  the  inviscid  veloc¬ 
ity  Induced  by  the  model,  thero  was  still  u  marked 


increase  in  the  magnitude  of  the  amplification  ae 
aora  blowing  wai  applied.  Accordingly,  the  present 
boundary- layer  computational  method  was  used  to 
examine  the  relationship  between  44  end  vw  ;  i.e., 
to  cast  light  on  the  amplification  of  the  normal 
velocity  acroe*  the  boundary  layer.  For  these 
calculations,  flow  without  tha  model  present  was 
represented  by  assuming  uniform  flow  over  a  flat 
surface  with  zero  straamwise  pressure  gradient. 

With  the  assumption  of  constant  transpiration 
velocity  along  the  controlled  section,  calculations 
wero  made  to  cover  the  range 

-0.0275  (  VM/Um  «  0.0150.  A  ot  of  V4/Uw  *14/4, 
over  part  of  this  range  is  given  Fig.  IS  for 
locations  corresponding  to  the  upstream  and  down¬ 
stream  limits  of  tho  controlled  test  section  length 
and  shows  only  a  weak  dependence  on  x  ,  When 
vw  /  U„  4  -0.0050,  thero  appears  t.o  be  an  asympto¬ 
tic  suction  behavior  of  the  same  nature  that  exists 
in  laminar  flow.35  That  is,  tho  boundary  luyer 
ceases  to  grow,  and  its  integral  properties,  parti¬ 
cularly  <5*  ,  become  independent  of  *  .  Eq.  (4) 
then  predicts  that  V*  •  vw,  so  thut  tho  slope  of 
tho  curvo  in  this  region  is  just  unity,  and  there 
is  no  amplification  by  the  boundary  layer.  Ilow- 
ovor,  as  V„/H,  increases  ubovo  -0.0050,  the  slope 
of  the  curve  in  Fig.  15  increasos  significantly, 
and  in  the  region  of  blowing,  v.  greatly  oxceods 
V^j  .  Not  only  has  the  First  term  on  the  right- 
hand  side  of  Uq.  (4)  Increased,  hut  this  in  turn 
has  groatly  magnified  tho  boundary- layer  growth 
represented  by  the  soci.nl  torm;  since  these  terms 
are  addltlvo,  umpl  1  f le.uti\  11  of  the  normal  velocity 
across  tho  layor  rosults,  .lacocks3(>  has  performed 
calculations  with  11  finitn-dlfforonce  boundary- 
layer  code  and  experimental ly-dotormlnod  wnll- 
pvossuro  dutu  with  n  model  present.  The  envelope 
which  he  dotorminod  for  his  calculated  results  is 
quite  similur  to  thut  exhibited  in  Fig.  15. 

In  uddition  to  helping  in  un  understanding  of 
tho  simulation  rosults  discussed  in  tho  main  text, 
Fig.  15  carries  a  very  practical  lesson  for  the 
users  of  such  adaptive-wall  tunnels;  vis.,  control 
of  tho  inviscid  flow  via  plenum  adjustments  is 
linear  und  well-behaved  for  suction,  but  likely  to 
bo  highly  nonlinear  when  blowing  is  uoedod. 


Append  I  x  11 

Idoa  1  lzod  Influence  Runet  I  ops 

lr  this  appendix,  an  Idealized  form  of  tho 
influence  functions  of  Section  IV  is  found.  A 
linear  approximation  for  the  flow  In  a  tunnel  with¬ 
out  thu  model  present  has  been  chosen.  Thus,  the 
idoalizod  analysis  boro  will  be  rolutod  to  the 
differences  found  in  tho  results  of  tho  simulations 
of  Suction  IV  when  the  pressure  Is  changed  in  ono 
plenum,  while  holding  tho  pressure  constant  in  tho 
other  chambers.  These  relationships  will  he  point¬ 
ed  out  as  tho  anulysis  proceeds. 

The  (  coordinate  system  used  here  is  tho 
same  us  the  X-  if  system  in  Fig.  1  except  that 
the  origin  of  %  is  located  in  lino  with  tho  center 
of  the  plenum  chamber  being  adjusted,  which  hus 
length  ,  whereas  the  origin  of  x  Is  at  tho 
quarter  chord  of  the  model.  In  tho  £  -  *1  systom, 
the  flow  investigated  is  induced  by  u  uniform 
normal  velocity,  ,  imposed  at  thu  upper  wull 
over  the  interval  -Xn/t  4  with  V  *  0 

prescribed  everywhere  else  on  the  upper  wull  us 
well  as  on  tho  ontire  lowor  wall.  The  froo-stream 


i) 


velocity  la  U„  and  the  etreamwlee  diaturbancaa , 
u.  t  are  assumed  to  be  zero  far  upstream,  juat  ;ia 
in  the  experiments.  The  magnitude  of  v,  ia  relat 
ed  to  the  differences  in  the  intagral  of  v«  (  *  ) 
over  the  length  of  planum  n  aa  carried  out  before 
and  after  changing  the  pressure  in  that  plenum. 


The  linearized  compreasible-flow  analysis  is 
carried  out  in  the  Prandtl-Glauert  approximation, 
so  that  a  disturbance  velocity  potential,  $  , 

exists  and  satisfies 

+  .■  0  (B-l 

where  fi* *  I  -  M*  and  the  subscripts  f  and  1 
denote  partial  differentiation  with  respect  to 
those  variables.  The  boundary  conditions  that  the 
u.  disturbance-velocity  component  must  be  zero  at 
infinity  upstream  and  that  v  at  the  walls  must 
be  as  described  above  are  expressed  at 
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Tho  problem  is  solved  by  meuns  of  Courier 
transforms,  which  aro  definod  as 
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4>  (*l  is;  «  ~~  J<p({,r\)e  d(  (B.S) 
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Upon  transformation,  Bq.  (B-l)  becomes 

5)  -  c  (b-7) 

which  in  an  ordinary  differential  equation  with  the 
solution 

/S  sij  -/}»n 

<p(n‘>v  *  +  fig  *  ft)  (b-8) 


Performing  the  appropriate  operations  on  Gqs.  (B-8), 
(B-9)  and  (B-10)  yields 

•-  „  *  V.  + 

1,9  "  (B.U) 
Ths  potential  itself  is  then  found  by  substi¬ 
tuting  in  Bq.  (B-6).  Since  the  right-hand  side 
of  Bq.  (B-U)  is  an  tvan  function  of  b  ,  this 
result  can  be  simplified  to 
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Differentiating  this  twice  with  respoct  to  £ 
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gives 
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which  can  bo  integrated^  to  give 


efih 

i 

CO^tv  | 

L  3/ih  \ 

4  Oi5  | 

TTCq+W 

.  zto  . 

AWvi-  ^ 

afi  h 

L 

j  +  0*4 

1  1 - 1 

-C 

* 

* 

(B-14) 


Tho  doslrcd  U  component  In  found  by  integra¬ 
ting  l!q.  (B-14)  with  respect  to  f  ,  immoly 

V*’#  “  ■  Trtib 


f  rcCf  +-t„/i)' 

+  cos 

r  ittythy 
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L  a  fih 

whore  is  ti  constant  of  integration.  This  Is 
evaluated  by  satisfying  tho  boundary  condition  in 
liq,  (B-2) ,  which  gives 


where  /?,  to  fl,  are  constants.  Sines  the  constant 
A|  transforms  back  to  the  coordinates  as  a 
delta  function,  which  has  no  physical  significance 
in  this  problem,  A,  is  taken  to  bo  zero.  The 
boundary  condition  of  Bq.  (B-2)  will  be  satisfied 
later  after  transformation  back  to  the  coordi¬ 
nate.  The  boundary  conditions  of  Bqs.  (B-3)  and 
(11-4 )  transform  to 

j  <- h ;a)  -  0  (B-9) 

^n(H;s)  *  2V,  »«(iIn/i)/VF?r  S  (B_i0) 
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so  that 
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The  desired  ir  component  is  foimJ  in  a  simi¬ 
lar  fashion.  Using  Bq.  (B-14)  in  conjunction  with 
Bq.  (B-i)  Rives  an  equation  for  (  df ,  yj  )  which 
can  be  integrated  with  respect  to  q  to  yield 
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where  the  conatant  of  integration  1»  »»  «*  *«• 
(B.S)  end  tjwUpis  *1  according  to  whether  f  *<»  . 

The  result*  are  now  cast  in  terms  of  the 
variables  *  and  y  .where  %  *  </-»'"  *"d 

j.n/N  .  In  particular,  the  idealised 

influence  functions,  denoted  by  1  *  *w!1*°  ***" 

tinguish  them  from  those  calculated  in  the  simula¬ 
tions,  are  found  from  Bqs.  (Bil7)  and  (B-18). 

From  Eq.  (B-17),  it  follows  that  ^ C -*•.#» J  *  1 0  *"d 
i,f<*?(|)«-v,Veh,so  that  in  terms  of  Bqs.  (5)  and 

(67.  S:*  •  Thu.  both 

and  must  be  divided  by  this  factor  to  get  1^, 
and  I  .  which  for  incompressible  flow  (/*•>)  are 
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Figs.  8%o  13  were  evaluated  from  BJj*‘  (B-19)T?"d 
CB-201  for  J »/h«  I  and  ii  *  *ic  »  t  a/S  .  The 
expression  for  iu«.V>  In  Eq.  (B-l'ilisvsUd 
everywhere  except  at  the  two  points  9  '»  *,  ’ 

where  there  are  singularities  introduced  by  the 
discontinuities  in  v  at  these  points.  A*  is  well 
known  in  linearized  thin-airfoll  theory,  such  a 
discontinuity  leads  to  a  logarithmic  singularity  in 
U  there.  (See  Refs.  38  and  39  for  the  singula¬ 
rity  at  the  kink  in  a  bent  tlat  plate.)  In  the 
calculated  values  of  l*.  shown  in  Figs.  B  and  10, 
there  are  no  singularities  because  i  *  */ '  ■  g* 

the  overshoot  observed  near  %  -t  ’/t  has  its  orig  n 
in  that  singular  behavior. 


•4.S7  -3.73  -3.7e"  .1.s«  -0M/ 


1.10  3.73  4.S7 

»/o 


,  TBAilUNRATKW  VILOCITY  AT  UPPER  TUNNEL  WALL  FOR 
,I,U"  '  MOMENTA FUINUM CONTROL.  NACA  0011  AIRFOIL, 

-  0.  ff  -  A  ,  V/®  ‘  10. 


11 


- UNCONPINKD  PLOW 

-  SMMINTSDPLSNUM  CONTROL 

- —  Acp.*o.a*3 

A  «  *  _  ft  M« 


Plturt  4  IPPICTIVI IN  VISCID  NORMAL  VILOCITY  AT  UPMR 

TUNNKL  WALL  DUS  TO  BOUNDARYLAYIR  DISP  LACS  MINT 
EPFICT  FOR  SEOMINTED-PLENUM  CONTROL.  NACA Mil 
AinPOtL, Mm  -0.  <f  "4  ,  v/«“1.0. 


__  _ 

_  — 

-  UNCONPINIO  PLOW 

X 

X 

X  8IOMINTIO  FLINUM  CONTROL 

Q 

0 

0  Acd  -°.»3 

O' 


•4J7  >3.71  .3.79  -1.M  0.11 


i-c.io 

u/U 

oo 

J^O.OI 

I-  0,04 
» 

r\ 

"°-02X» 

.... 

0,1?  1.47  2.19  3.73  4.97 


«/« 


UPPSR  CONTROL  iUNPACI.  y ./«  •  0.8M 


1.39  2,39  3.94  4.17 


1—9,99 

LOWER  CONTROL  SURPACI,  v,/«  ■  0.494 

Plfuti  S  COMPARISON  OP  STRSAMWISI  OISTURBANCI  VILOCITY  AT 
CONTROL  SURPACIS  POR  UNCONPIKID  PLOW  AND 
I1MMNTBD-PLBMUM  CONTROL.  NACA  Mil  AIRPCHL, 

Mm  «0,  <T  "*  4*. 


- UNCONPINIO  PLOW 

X  X  X  SSOMSNTID-PLSNUM  CONTROL 
o  O  O  top  -4AM 


LOWIN  CONTROL  SURPACI,  v,/a  •  4.4(4 

PifUrt  S  COMPARISON  OP  NORMAL  VELOCITY  AT  CONTROL 

SURFACES  FOR  UNCONFINED  FLOW  ANOIEQMENTED- 
PLENUM  CONTROL.  NACA  0012  AIRFOIL,  »„  -0,  rr  .  4° 


2.0 


r 
.1.81*1 


- - - UNCON  PI  NED  FLOW 

X  X  X  «OM£NTBD-rUNUM  CONTROL 


i.oh  x 


O.Dh 


"7 


1.0 


\ 


"‘X-. 


’ '  *-K*. , 


..-,,Li.d;.r,)l''r:>l  v~-K' -in  ■«■■■. — — i 


0.2 


0.4  0,4  0.4"*';-'ft.-jL9 

k/« 


Flour*  7  COMPARISON  OP  AIRPOI L  FRESSURS  DISTRIBUTIONS  FOR 
UMCONFINID  PLOW  AND  IIQMINTIOPLENUM  CONTROL. 
NACA  Mil  AIRFOIL,  M,„  »0.  nr  -  4°, 


12 


—  IDEALIZED?,.  ,  APPENDIX  I 

X  NUMERICAL  BHMJLATIOMt 
WITH  MODEL 


- 


,  BOUNDARIES  Of  PLENUM 
I  SEINE  ADJUSTED 


•OUNDANIII  Of  PLENUM  S  / 


M 


FIbuk  •  INC LUINCE  FUNCTION!  FOB  ITNCAMWIU  VELOCITY 
COMPONENT  AT  CONTROL  fUNFACI  ADJACENT  TO 
PLENUM  BEINQ  ADJUSTED, 


Flour.  10  AMPLIFIED  DIVTRIBUTIONI  OF  INFLUENCE  FUNCTION!  FOR 
ITRIAMWIIE  VELOCITY  COMPONENT  AT  CONTROL  IURFACE 
ADJACENT  TO  PLENUM  A 


- IDEALIZED^  .  APPEND!  A  I 

I 

X  X  NUMERICAL  SIMULATIONS 
WITH  MODEL 


IDEALIZED*,  ,  APPENDIX  I 


BOUNDARIES  OP  PLENUM 
■SINO  ADJUSTED 


-o —  -0.0E3 

u 

--■a —  Ac„  --0.031 
**• 


BOUNDARIES  OP  PLENUM  E 


ar**  — 

f 

!  0.1 

Fli  ir.  •  INFLUENCE  FUNCTION!  FOR  NORMAL  VELOCITY 
COMPONENT  AT  CONTROL  IURFACE  ADJACENT  TO 
PLENUM  HIND  ADJUITED. 


Flsurt  11  AMPLIFIED  DISTRIBUTION!  OF  INFLUENCE  FUNCTION! 
FOR  NORMAL  VI  LOCI  TY  COMPONENT  AT  CONTROL 
IURFACE  ADJACENT  TO  PLENUM  A 


DISTRIBUTION  LIST  TOR  URCI/-S8ITIZD 
TECHNICAL  REPORTS  AND  REPRINTS  ISSUED  UNDER 

conotact  t/OCMf-WrA  -jfttfmsK 

All  addresses  receive  on*  copy  uni***  oth«rvl««  specified 


Technical  Library 
Building  313 

Ballistic  R****rob  Laboratories 
Aberdeen  Proving  Ground*  MD  21005 

Mr.  Aviar*  Celains 
Balliatic  R***arch  Laboratory 
Ballistic  Modelling  Division 
Ab*rd**n  Proving  Ground*  MD  21005 

Dr.  P.  J.  Roach* 

Ecodynemia*  Research  Associates*  Inc. 

P.  0.  Box  8172 
Albuq.uerq.ue  *  KM  8T106 

Defense  Technical  Information  Center 

Cameron  Station,  Building  5 

Alexandria,  VA  22314  12  copies 

Library 
Naval  Academy 
Annapolis,  MD  21402 

Director,  Tactical  Technology  Office 
Defense  Advanced  Research  Projects 
Agency 

1400  Million  Boulevard 
Arlington,  VA  22209 

Cod*  211 

Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington*  VA  22217 

Cod*  438 

Office  of  Naval  Researoh 
800  N.  Quincy  Streat 

Arlington,  VA  22217  2  conies 

Dr.  J.  L.  Potter 
Deputy  Director,  Technology 
von  Kerman  Gas  Dynamics  Facility 
Arnold  Air  Force  Station,  TN  37389 

Profess  r  J,  C.  Wu 
School  '  f  Aerospace  Engineering 
Georgia  Institute  of  Technology 
Atlanta  QA  303,32 


Library 

Aerojet-Genera1  Corporation 
6352  North  Irv*^iale  Avenue 
Azusa*  CA  91702 

NASA  Scientific  and  Technical 
Information  Facility 
P.  0.  Box  8757 

Baltlaore/Vashington  International 
Airport »  MD  21240 

Dr.  K.  C.  Wang 
Martin  Marietta  Corporation 
Martin  Marietta  Laboratories 
1450  8outh  Rolling  Road 
Baltimore*  MD  2122T 

Professor  A.  J.  Chorln 
Department  Mathematics 
University  of  California 
Berkeley,  CA  94720 

Professor  M.  Holt 

Department  of  Mechanical  Engineering 
University  of  California ’ 

Berkeley,  CA  94720 

Dr.  H.  R.  Chaplin 
Code  1600 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Betheeda*  MD  20084 

Dr.  Hans  Lugt 
Cod*  l602 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Betheeda*  MD  20084 

Dr.  Francois  Frtnkiel 
Cods  1602 

David  W.  Tsylor  Naval  Ship  Research 
and  Development  Center 
Betheeda,  MD  20084 


Pams  2 


Dr.  0.  R.  Inger 

Department  of  Aerospace  Engineering 
Virginia  Polytechnic  Institute  and 
State  University 
Blacksburg,  VA  2U061 

Professor  C.  H.  Levis 
Department  of  Aerospace  and  Ocean 
Engineering 

Virginia  Polytechnic  Institute  and 
State  University 
Blacksburg,  VA  2U06l 

Professor  A.  H.  Hayfeh 
Departnent  of  Engineering  Science 
Virginia  Polytechnic  Institute  and 
State  university 
Blacksburg,  VA  24061 

Dr.  A.  Rubel 
Resevroh  Department 
Grumman  Aerospace  Corporation 
Bethpage,  HI  11714 

Commanding  Officer 
Office  of  Naval  Research  Eastern/ 
Central  Regional  Office  (Boston) 

666  Summer  Street,  Bldg.  114,  Section  D 
Boston,  MA  02210 

Dr.  J.  C.  Erickson,  Jr. 

CAL8PAH  Corporation 
Advanced  Technology  Center 
P.  0.  Box  400 
Buffalo,  HZ  14225 


Mr.  C.  E.  Wittliff 
CALSFAH  Corporation 
Advanced  Technology  Center 
P.  0.  Box  400 
Buffalo,  HZ  1422$ 

Professor  R.  F.  Probstein 
Department  of  Mechanical  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge,  Ml  02139 


Commanding  Officer 

Si?1®*  of  Have!  Research  Branch  Office 


South. 


PUS? 


Coda  753 

Haval  Weapons  Csnter 
China  Lake,  CA  93555 

Mr.  J.  Marshall 
Code  4063 

Haval  Weapons  Center 
China  Lake,  CA  93555 

Professor  R.  T.  Davis 
Department  of  Aerospace  Engineering 
University  of  Cincinnati 
Cincinnati,  OH  45221 

Professor  S.  0.  Rubin 
Department  of  Aerospace  Engineering 
sad  Applied  Mechanics 
University  of  Cincinnati 
Cincinnati,  OH  45221 

Library  MB  60-3 
RASA  Levis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 

Dr.  J.  D.  Anderson,  Jr. 

Chairman,  Department* of  Aerospace 
Engineering 
College  of  Engineering 
University  of  Maryland 
College  Park  MD  207U2 

Professor  0.  Burggraf 
Department  of  Aeronautical  and 
Astronautical  Engineering 
Ohio  State  University 
1314  Kiaaear  Road 
Columbus ,  OH  43212 

Technical  Library 
Haval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren,  VA  22446 

Dr.  F.  Moore 

Haval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren ,  VA  22446 

Technical  Library  2-51131 
LTV  Aerospace  Corporation 
P.  0.  Box  5907 
Dallas,  TX  75222 


fjfiT.  3 


Library*  Unit ad  Aircraft  Corporation 
Research  Laboratories 
Silver  Lane 

Kaat  Hartford,  CT  06108 

Professor  0.  Noretti 
Polytechnic  Inatltuta  of  Hew  York  Lone 
Island  Center 

Department  of  Aerospace  Engineering  and 
Applied  Mechanics 
Route  110 

framing dale,  BY  11735 
Dr.  V.  E.  Briley 

Scientific  Research  Associates,  Inc. 

P.  0.  Box  1*98 
Glastonbury,  CT  06033 

Professor  P.  Gordon 
Calumet  Campus 
Department  of  Mathematics 
Purdue  University 
Haasond,  IB  1+6323 

Library  (MS  183 ) 

HASA  Langley  Research  Center 
Langley  Station 
Hampton,  VA  2366s 

Professor  A.  Chapsumn 
Chairman,  Mechanical  Engineering 
Department 

Villiam  M.  Rice  Institute 
Box  1892 

Houston,  TX  77001 

Technical  Library 
Haval  Ordnance  Station 
Indian  Head,  MD  2061+0 

Professor  D.  A.  Caughey 
Sibley  School  of  Mechanical  and 
Aerospace  Engineering 
Cornell  University 
Ithaca,  HY  14850 

Professor  E.  L.  Resler 
Sibley  School  of  Mechanical  and 
Aeroapace  Engineering 
Cornell  University 
Ithaca,  NY  14850 


Professor  8,  f,  Shen 
Sibley  8chool  of  Mechanical  and 
Aerospaoe  Engineering 
Cornell  University 
Ithaca,  BY  14850 

Library 

Midwest  Research  Institute 
425  Volkar  Boulevard 
Kansas  City,  M0  64uo 

Dr.  M.  M.  Safes 
Flow  Research,  Inc. 

P.  0.  Box  5040 
Kent,  VA  98031 

Dr.  E.  M.  Muraan 
Flow  Research,  Inc. 

?.  0.  Box  5040 
Kant ,  VA  98031 

Dr.  J.  J.  Riley 
Plow  Research,  Inc. 

F.  0.  Box  5040 
Kent,  VA  98031 

Dr.  S.  A.  Orszag  . 

Cambridge  Hydrodynamics,  Inc. 

54  Baskin  Road 
Lexington,  MA  02173 

Dr.  P.  Bradshaw 

Imperial  College  of  Science  and 
Technology 

Department  of  Aeronautics 
Prince  Consort  Road 
London  SV7  2BY ,  England 

Professor  T.  Cebeci 
Mechanical  Engineering  Department 
California  State  University,  Long 
Baach 

Long  Beach  ,  CA  90840 
Dr.  H.  K.  Cheng 

University  of  Southarn  California, 
Dapartmant  of  Atrospaca  Engineering 
University  Park 
Los  Angeles,  CA  90007 


Professor  J.  D.  Cols 
Mechanics  sad  Structures  Department 
School  of  Engineering  sad  Applied 
Science 

University  of  California. 

Los  Angeles (  CA  90024 

Engineering  Library 
University  of  Southern  California 
Box  T7929 

Los  Angeles,  CA  90007 
Dr.  C.  -M.  Ho 

Department  0f  Aerospace  Engineering 
University  of  Southern  California 
University  Park 
Los  Angeles,  CA  90007 

Dr.  T.  D.  Taylor 

The  Aerospace  Corporation 

P.  0.  Box  92957 

Los  Angeles,  CA  90009 

Cosuoandlng  Officer 
Naval  Ordnance  Station 
Louisville,  KY  4 021 4 

Mr.  B.  H.  Little,  Jr. 
Lockheed-Qeorgla  Company 
Department  72-74,  Zone  369 
Marietta,  QA  30061 

Professor  E.  R.  Q.  Eckert 
University  of  Minnesota 
241  Mechanical  Engineering  Building 
Minneapolis ,  MH  55455 

Dr.  Gary  Chapman 
Mail  Stop  227-4 
Ames  Research  Center 
Moffett  Field,  CA  94035 

Library 

Naval  Postgraduate  School 
Monterey,  CA  93940 

Dr.  J.  L.  Stager 
Flow  Simulations,  Inc. 

735  Alice  Avenue 
Mountain  Viev,  CA  94041 

Dr,  S.  8.  Stabara 

Uijl*en  Engineering  and  Research 

£2n?iiS\ftr!A 

94043 


Engineering  Societies  Library 
345  East  4Tth  Street 
New  York,  NY  10017 

Professor  A.  Jameson 
C  our  ant  Institute  of  Mathematical 
Sciences 

New  York  University 
251  Mercer  Street 
New  York,  NY  10012 

Profee sor  G.  Miller 
Department  of  Applied  Soienee 
Nav  York  University 
26-36  Stuyvesant  Street 
New  York,  NY  10003 

Office  of  Naval  Researoh 
Nav  York  Area  Office 
715  Broadvay  -  5th  Floor 
New  York,  NY  10003 

Dr.  ft  Vagllo-Laurin 
Department  of  Applied  Science 
New  York  University 
26-36  Stuyvessnt  Street 
Nev  York,  NY  10003 

Mr.  D.  Farmer 

Naval  Ocean  Research  and  Development 
Activity 
Code  332 

N8TL  Station,  MS  39522 

Librarian,  Aeronautical  Library 
National  Research  Counoil 
Montreal  Road 
Ottawa  7,  Canada 

Lockheed  Missiles  and  Space  Company 
Technical  Information  Center 
3251  Hanover  Street 
Pale  Alto,  CA  94304 

Commanding  Officer 
Office  of  Naval  Research  Western 
Regional  Office  (Pasadena) 

1030  Bast  Green  Street 
Pasadena,  CA  91106 

Engineering  Division 
California  Institute  of  Technology 
Pasadena,  CA  91109 


Pam  5 


Library 

Jat  Propuls ioa  Laboratory 
1*800  Oak  Or  ova  Drivt 
Pasadena ,  CA  91103 

Professor  H.  Lispmann 
Department  of  Aeronautics 
California  Instituts  of  Technology 
Pasadana ,  CA  91109 

Mr.  L.  I.  Chaaau,  MQR-MBD  Lib. 

Oanaral  Electric  Company 
Misalla  and  Spaoa  Division 
P.  0.  Box  8555 
Philadelphia,  PA  19101 

Tachnioal  Library 
Naval  Missila  Cantar 
Point  Mugu,  CA  93042 

Profassor  S.  Bogdonoff 
Oas  Dynamics  Laboratory 
Department  of  Aarospaca  and 
Mechanical  Sciancas 
Prinoaton  Univarsity 
Princeton i  NJ  08540 

Professor  8.  I.  Chang 
Departmant  of  Aerospace  and 
Maohanioal  Sciancas 
Prinoaton  University 
Princeton,  NJ  08540 

Dr.  J.  E<  Yates 

Aeronautical  Research  Associates 
of  Princeton,  Inc. 

50  Washington  Road 
Princeton,  NJ  08540 

Profassor  L.  Sirovich 
Division  of  Applied  Mathematics 
Brown  University 
Providence ,  RI  02912 

Redstone  Scientific  Information  Center 
Chief,  Document  Section 
Army  Missile  Command 
Redstone  Arsenal,  AL  35809 

U.S.  Army  Research  Office 

P.  0.  Box  12211 

Research  Triangle ,  NC  27709 


Editor,  Applied  Mechanics  Review 
Southwest  Re search  Institute 
8500  Culebra  Road 
San  Antonio,  TX  78226 

Library  and  Information  Services 
Otneral  Dynamics -CO NVAIR 
P.  0.  Box  1128 
San  Disgo,  CA  92112 

Dr .  R .  Magnus 
General  Dynamios-CONVAXR 
Kearny  Mesa  Plant 
P.  0.  Box  80647 
San  Disgo,  CA  92138 

Offioa  of  Naval  Research 
San  Francisco  Area  Offioe 
One  Hallidie  Plaza,  Suite  601 
San  Francisco,  CA  94102 

Library 

The  RAND  Corporation 
1700  Main  Street 
Santa  Monica,  CA  90401 

Dr.  P.  E.  Rubbert 
Boeing  Aerospace  Company 
Boeing  Military  Airplane  Development 
Organization 
P.  0.  Box  3707 
Seattle,  WA  96124 

Dr.  K.  Yoshibara 
Boeing  Aerospace  Company 
P.  0.  Box  3999 
Mail  Stop  41-16 
Seattle ,  WA  96124 

Librarian 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 

Dr.  J.  M.  Solomon 
Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  2<*910 

Professor  J.  H.  Ferzlger 
Departmant  of  Mechanical  Engineering 
Stanford  University 
Stanford,  CA  94305 


?mbJl 

Professor  K.  Karwacheti 
Department  of  Aeronaut let  and 
Astronautics 
Stanford  university 
Stanford,  CA  94305 

Professor  0,  Bunenana 
Institute  for  Plasna  Research 
Stanford  University 
Stanford,  CA  94305 

Engineering  Library 
McConnell  Douglas  Corporation 
Department  218,  Building  101 
P.  0.  Bo*  51 6 
St.  Louis,  M0  63166 

Dr.  R.  J.  Hakkinen 

McDonnell  Douglas  Corporation 

Department  222 

P.  0.  Bo*  516 

St.  Louis,  MO  63166 

Dr.  N.  Malmuth 

Rockwell  International  Science  Center 
1049  Camino  Dos  Rios 
P.  0.  Bo*  1085 
Thousand  Oaks,  CA  91360 

Library 

Institute  of  Aerospace  Studies 
University  of  Toronto 
Toronto  5,  Canada 

Professor  W.  R.  Sears 

Aerospace  and  Mechanical  Engineering 

University  of  Arizona 

Tucson,  AZ  85721 

Professor  A.  R.  Setbass 
Department  of  Aerospace  and 
Mechanical  Engineering 
University  of  Arizona 
Tucson ,  AZ  85721 

Dr.  K.  T.  Yen 
Code  3015 

Naval  Air  Development  Center 
'Werminater,  PA  18974 


Air  Force  Office  of  Scientific 
Research  (SR3N) 

Building  l4l0,  Bolling  AFB 
Washington,  DC  20332 

Chief  of  Research  and  Development 
Office  of  Chief  of  Staff 
Department  of  the  Army 
Washington ,  DC  20310 

Library  of  Congress 

Scienoe  end  Technology  Division 

Washington,  DC  20540 

Diractor  of  Research  (Code  RH) 
National  Aeronautics  and  Space 
Administration 
600  Independence  Avenue,  SW 
Washington,  DC  20546 

Library 

National  Bureau  of  Standards 
Washington,  DC  20234 

National  Scienos  Foundation 
Engineering  Division 
3.800  C  Street,  NW 
Washington,  DC  20550 

AIR  32CD 

Naval  Air  Systems  Command 
Washington,  DC  20361 


AIR  950D 

Naval  Air  Systems  Command 
Washington,  DC  20375 

Cods  2627 

Naval  Research  Laboratory 
Washington,  DC  20375 

SEA  03512 

Naval  Sea  Systems  Command 
Washington,  DC  20362 

SEA  0903 

Naval  Sea  Systems  Command 
Washington,  DC  20362 


». 

L 


gMMLl 


Dr.  A.  L.  Slafkosky 
Scientific  Advisor 
Cowandant  of  the  Marin*  Corps 
(Cod*  AX) 

Washington,  DC  20380 
Director 

Weapons  Systems  Evaluation  Group 
Washington,  DC  20 350 

Research  Library 
AVCO  Corporation 
Missile  Systems  Division 
201  Lovsll  Street 
Wilmington,  MA  0168? 

AFAPL  (APRC) 

AB 

Wright  Patterson,  AFB,  OH  1*5433. 

Dr.  Donald  J.  Harney 
APPDL/PX 

Wright  Patterson  APB,  OH  1*51*33 


Dr.  D.  R.  Kotansky 
Department  34  J, 

P.  0.  Box  S16 

St.  Louis,  M0  63166 


